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NASA Grant Number NGR 34-002-193 
SURFACE FITTING THREE-DIMENSIONAL BODIES 

by Fred R. DeJarnette 
North Carolina State University 

SUMMARY 

The geometry of general three-dimensional bodies is generated from coordinates 
of points in several cross sections. Since these points may not be smooth, they 
are divided into segments and general conic sections are curve fit in a least-squares 
sense to each segment of a cross section. The conic sections are then blended in 
the longitudinal . direction by fitting parametric cubic-spline curves through coordi- . 
nate points which define the conic . sections in the cross-sectional planes o Both 
the cross-sectional and longitudinal curves may be modified by specifying particular 
segments as straight lines and slopes at selected points. Slopes may be continuous 
or discontinuous and finite or infinite. 

After a satisfactory surface fit has been obtained, cards may be punched with 
the data necessary to form a geometry subroutine package for use in other computer 
programs. At any position on the body, coordinates, slopes and second partial 
derivatives are calculated® 

ic til 

The method is applied to a blunted 70 delta wing ^ and it was found to generate 

eis.Lmo 

the geometry very wellc 



INTRODUCTION 


Many disciplines require a mathematical description ei 
surfaces which cannot be represented by simple mathematical expressions r. The 
location and slopes of points on a body are needed in the analysis or jjcriictct ee 
and inviscid flow fields^ Viscous flow-field analyses require the body cur/atuoa 
in addition to location and slopes o The geometrical properties of simple shapes 
like spheres > cones, ellipsoids ^ and paraboloids can be described by relatively 
simple mathematical equations o However, many configurations of interest today, 
such as the space shuttle, are complex three-dimensional shapes whose geometry 
cannot be described easily., In mahy instances all the information that is given 
is a drawing with a plan view, a side view, and several cross sections of the 
vehicle^ Sometimes models of a vehicle are available, and coordinate positions 
can be accurately measured on them^ On the other hand, slopes cannot be accurately 
measured on a drawing or a model, and the determination of radii of curvature is 
even less accurateo This report develops a computer program which will surface 
fit mathematical relations to complex three-dimensional bodies , The method 
yields accurate coordinates and slopes and reasonably accurate radii of curvature 
at any position on the body„ 

Previous approaches to surface fitting three-dimensional bodies generally 
divided the surface into **patches*' and represented each patch by flat surfaces, 
cufelc or higher order polynomials, or conic sections (eogo, see refs, 1-3), 

Each of these methods has undesirable features in the forms used previously. 

Flat surfaces are completely unacceptable if radii of curvature or continuous 
slopes are needed. Cubic or higher order polynomials often lead to unwanted 
wiggles and bulges since they allow points of inflections. In ref, 3 this author 
applied the method of double splines, which used bi-cubic interpolation, to 
surface fit the coordinates of points in several cross-sectional planes of 
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three-dimensional bodies « However, It was. found that bulges and/or dimples 
occurred in these surfaces, partlcu],arly when the thickness was_ much smaller 
than the. spane ' 

Reference 4 approximated the shape of a space shuttle orbiter with elliptical, 
cross sections with different elllpticity on the windward and leeward sideso 
Cubic polynomials were used to define segments of the plan and thickness distribu- 
tions c The coefficients of these polynomials were .chosen to make the slopes 
continuous across boundaries of the longitudinal segmentSo However, points of 
inflections were found to occur inside the segments- and thus gave the undesirable 
bumps, which in turn significantly affected the surface pressure and heating 
distributions calculated with this geometry,, . 

Cpnio sections have been used to describe segments of aircraft contours for 
some time (see ref o ..1) » Both longitudinal and transverse contours were represented 
by conic ' sections , but the slopes at the ends , of each segment had to be measured » 

As mentioned earlier, slopes are difficult to measure accurately. The present . 
author developed a least-square curve-fitting technique (ref » 5) using general 
conic sections for each segment .of a curve in a plane.,, .Since all data points are 
not necessarily smooth, the curve was constrained to go through selected control 
points (which are boundaries for each conic . aegmept) but fit the remaining data 
only in a least-square sense. In addition,, the .slope at each control point was 
constrained to be continuous. The results from this curve-fitting method were , 
quite good, ' Another feature of conic sections is that they cannot yield curves 

J,. . ‘ . 

i 

with Inflectioh'points^ and this feature can.be a definite advantage. If an, 
inflection point .should be needed, then it could be made a control point with the 
resulting curves- free of inflection points except at this control .point, • 

Coons (ref, < 2) developed a sophisticated technique to describe three-dimensional 
surfaces by using blending functions to blend the surf ace between the boundary 
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curves of each patcho A major difficulty in applying Coon's method is that 
the user must supply the coordinates » slopes « and twists (cross derivatives) at 
all four corners of each surface patch* This information is generally difficult 
to determine, particularly the cross derivatives (see ref* 6)» Another disadvantage 
of using Coons’ method in computer applications is that 64 parameters are needed 
to describe each surface patch, and for a large number of patches the amount of 
storage required by the computer may be excessive* 

The method presented here for surface fitting three-dimensional bodies is 
somewhat similar to Coons’ patching method, but it reduces the number of parameters 
required to. describe a surface patch and is simpler to apply* Data points in 
cross sectional planes are curve-fit in a least-squares sense by segments of general 
conic sections* The conic sections are then blended in the longitudinal direction 
by fitting parametric splines through coordinate points which define the conic 
sections in the cross-sectional plapes* This technique has the advantage of 
allowing the user to continually modify the cross-sectional curves and the 
longitudinal curves until the body shape has the desired features* . Discontinuous 
slopes in both the circumferential and longitudinal directions may be specified. 
After the body shape has been modified to the final form, the variables 
necessary tq describe the shape may be punched on data cards*. Then for applications 
of the method, a relatively small geometry subroutine can be used along with these 
data cards to calculate positions, slopes, and radii of curvature on three- 
dimensional bodies* 
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paranecer defined by eq» (26) 

constant vector defined by eq« (20) 
body radius deflnedr, by,. eq» , (22) . 

r '• A. .rsi* i.' . 'V , 

matrix defined, by eq« (21) 
parameter defined by eq* (27) 

data point number of first control point in segment j 

slopes defined by eqsa (5) and (6) 

number of segments in a cross section 

terms defined by eqse (A2) and (A3) 

residual in eqs, (Bl) and (B2) 

term defined by eq» (AS) 

chordal distance bet\?een coordinates. 

parameter defined by eq« (25) 

Cartesian coordinates ^ see Figure 2 
local coordinates, see Figure 3 
coefficients defined by eqs« .(15), (16), and (17) 
slope of segment j , see Figure 5 





circumferential angle defined by eq» -(23) 


Subscripts 

h Intermediate point in a cross.^8ectlonal segment 
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j segment number In a cross section 

k data point number In a cross section 

r reference point 

s slope point of a cross-sectional segment 

0 first control point of a cross-sectional segment 

1 last control point of a cross-sectional segment 


Superscript 

* 


value where both R 


J 


0 and dRj / dy =0 
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ANALYSIS 

a tTpical application^ the geometry. of a body must be determined from a. 
model;, such as a vind tunnel models or a . three-view drawinge Generally several 
cross .sections , are obtainable from the. model or drawings the coordinates of 
data points on the boundary -of these cross 'Sections ' can be measured or calculated c 
A ^three-dimensional surface must then be fit through all the body cross sections n 
Many times the data points , in the cross sections are not completely smooths 
in those cases a smooth surface cannot be fit through all the data pointSe It 
is. then desirable to, have a smooth surface pass through designated data points ^ 
called control points » and pass, close tOp but not necessarily through, the other 
data points d v 

In thS: method presented heres- the data . points in each cross section are 
divided into segments, and portions of general ..conic sections are curve-fit to 
the data points .In each segment « ' The data points at ^he ends of each segment 
are designated as control points, and the curve ,1s constrained to go through. the 
control points as shown, in Figure It A three-dlmenslpnal surface le then generated 
by ’’blending” • the cross-sectional curves In the longitudinal direction c , Consider 
first the. technlciue for curve-fitting the data points In a cross-sectional planes 


Curve-Pit in a CrpsSTSectipnal Plane 
As mentioned previously, the; data points, in a crossr-sectional plane are 
generally npt completely smooth, and in. those cases a smooth curve cannot. be 
made to ,pasp .through every data .point. ’ Therefore, the data points are divided 
into segments, and a, portion of a general conic section is curve-fit in a leasts 
squares .sense through the data points in that segment,. The curve is constrained 
to go through the control points (end points, pf a segment) apd also have a 
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continuous slope at each control point unless the slope at a control point is 
specified otherwise o 

Define a three-dimensional coordinate system y» 5 with x in the longitudi- 
nal direction and x = constant is a cross-sectional plane (see Figure 2)<. Let j 
denote the, segment number in a cross-sectional plane as shown in Figure I 9 with 
the first segment starting on the positive y - axis and j Increasing clockwiseo 
For each segment, it is convenient to initially use a local coordinate system y^z 
with the origin on the first control point and the positive y-axis passing through 
the control point at the other end of the segment (see Figure 3)» In this local 
coordinate system.lt is easy to investigate possibilities of complex roots and 
interpret the coefficients of a general conic section geometrically » whereas it 
is difficult to interpret them geometrically in the global coordinates y^Ze 

For the j segment the equation for a general conic section is given by (refo7) 

A^y^ + B^yz + Cj,z^ + D^y + E^z + F^ =0 (1) 

Only 5 of the 6 coefficients are independent since the equation may be divided 
by any non-zero coefficient 6 The constraints that the curve pass through the two 
control points 

y “ 0 , z “ 0 and V “ Vj » z “ 0 yield 

F. » 0 (2) 

J 

Dj - - (3) 

By differentiating eq. (1), the elope In local coordinates is given by 


dz 

dy 






V 


- 2A,y - B^z 



+ 2C z + E. 

i J 


(4) 
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Although the slopes ,, at the ends of the segments are generally not known p define 




at y " 0 ^ z “ 0 


(S) 


and 



at y ■ yj t z “ 0 


( 6 ) 


Using these, two equations in eqo (4) it 'fol^.Sws that 


^3 ■ ^ ^3 '“3 

(?) 

and 


B «-A(^ + ^) 
j j mj 1^ 

00 

These results show that the end slopes do not affect the coefficient C^o On the 
other hand, the product A^Gj determines the, nature of the general conic .section 

A 1^ 

(refs -7)o If ^ 4A^ ® ,0 the conic is a parabola^ if -- 

4AjCj 0 the conic 

2 

Is an ellipse » and if « 4AjCj> 0 the conic. is a hyperbplao 

Of particular 


Interest la the poBslbility of complex roots when solving for z as a function of 
y in the region of interests It .is shown in Appendix A that z will have no 
complex -roots in the region 0 < y < y^ if 


AjCj..,-' 


For prescribed slopes m^. 


n^ it is interesting to observe ,, how the product AjCj 


affects a conic- section '(See Figure 4) » 


Unless a slope is specified at ^a control point, the method used here constrains 
the . slope in global coordinates :(y. z). to be continuous at a control point o This 

makes the. conic section in one segment dependent on data points in other segments 
as well as its own o From Figure 5 it can he seen that continuity of slope at 
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control point j requires 


tan”^ n^__^ “ tan~^ 


This equation can be expanded and learranged in the form 


ao> 


sin 46 ^ cos A6^ 


J“1 


(A/mJ^ cos 46 . - A. sin 40. 
3 3 j j 


Cil'^ 


This last equation is non-linear in the coefficients A^ » which would cause 
difficulty, in obtaining a solution for them. However $ as mentioned previously 
only 5 of 6 coefficients -in a given segment are independent c This allows an 
additional constraint to be imposed without affecting the overall equation for 
the general conic section. The additional . constraint used here is to equate 
separately the numerators and denominators of both sides of eq? (11). This gives 
two equations which are linear in the coefficients A^ ^ and which may be combined 
to. yield 


(A/n)j “ (AjCos40j^j^ "^j+1^ 


(12) 


(A/m) 


. “ (A. --A.COS4Q.) /sin40. 

J j”l J J i 


(13) 


For t;he first segment (j “ 1) m^ is specified and, equation (13) is not needed^ 
and also for the last segment (J » N) n^ is specified and eq. (12) is not needed; 

Now substitute eqs. (2), (3)> (7), (8), (12), and (13) into eq. (1) to 
obtain the equation for the conic sections as 




(14J 
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where 


(zyj-2y>/sinA6^ for j > 1 


0 for j =1 


y + cotA 02 )yz + y^z/jn^^-y^y for j ® 1 

2 

y +(cotA6j“CotA9j^^)yz - cotASjy^z-y^y for 1 < j < N 
2 

y +(cotA0jj”l/r^)yz - cotAe^^y^jZ - y^^y for j » N 


yz/sinA0^^^ 

0 for j ».N 


for J < N 


These results indicate that there are only two unknowns in each segment ^ Aj 
and However;, there. is one more segment than interior control points which 

means one of these coeffi^cletlta is arbitrary® Here » 1 is used arbitrarily 
(unless the conic section .requires A^^ = 0) and the independent unknowns become- 
» Aj (j N) o ■. If there were as many data points as nnknownsp and 

if some conic section could be made to pass through these data points ^ then 
eq®.(14) could be- applied to all the data points to give (2N"1) linear equations 
for the (2N”1) coefficients C^ip p A^ for j “ 2g'‘“o.p N® However p there are 
generally more than (2N'-il) data points* and if eq® (14) were applied to all of 
them an'overdetermined system of linear equations (see ref® 8) would resplt® 
Therefore* a . least-squares solution of the overdetermined system is used to 
determine the coefficients C^* Cj* Aj for j =» 2*,“'>% N, This procedure ,1s 
described in Appendix -B® After obtaining the coefficients in this manner p the 
inequality of eq® (9) is checked. for the possibility of complex roots for z in 
each segment® If the inequality of eq® (9) is not satisfied* the coefficient C.- 
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in that section is replaced by the value obtained using the equality sign in 
eqo (9 ) 0 As mentioned earlier^ the coefficient does not affect the slopes of 
the curve at the end points of the segment ^ and hence does not affect other 
segments^ The equality sign in eq» (9) gives 2 straight lines as shovn in 
Figure 4o 

In order to solve for z as a function of y from eq, (14) > a quadratic 
equation must be solved and the proper choice of the + or - sign must be 
determined beforehand for each segment. It is shown in Appendix A that in order 
to make z s 0 at y = y^ (a control point at the end of segment j)» the + sign 
must be used if (A/m). > 0 and (A/n). < 0 ^ and the - sign must be used if (A/m).<D 

J j 1 

A 

and (A/n)^ > Oi. 

Once the coefficients Aj and are determined, all the conic sections are 
completely defined for a given cross-sectional plane. In order to put these 
results iifito a form suitable for "blending" the cross sections in the longitudinal 
direction, the conic section for each segment is redefined in terms of 4 points 
- the two control points at the ends of the segment, a slope point which determines 
the slopes at the end points, and finally an intermediate point on the curve 
between the end points (see Figure 6). The 3 points on the curve and the two 
slopes at the end points of a segment are sufficient to determine new coefficients 
A^, A 2 , A^, A^ , and A^ for the general conic section, 


A^y^ + A^yz + A^z^ + A^y + A^z + 1 


(18) 


*Note that in order for a conic section to pass through the control points, 

n. <0 if m. > 0 and n. > 0 if m^, < 0. 

3 J j j 

**Note that the conic section given by eq« (18) must have the constant 1 replaced 
by 0 if the curve is to pass through the origin y « 0 , z 4 0. 
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in global coordinates y^z (see re£o 1) o This process is done for each segment 
in a cross section, and the 5 coefficients will become functions of the icngitudinal 
coordinate x_^hen blending a. segment pf_a„,cross eeetion .with corresponding segments 
in the other cross sections. 

Longitudinal Variation of Cross Section 

In order to determine the longitudinal variation of -the coefficients in 
eq, (18), a three-dimensional curve is fit through each of the 4 points used to 
define the conic section of corresponding segments (see Figure /)o In contrast 
to the cross sectional data points, these curves must pass through each ©f the 
points in the longitudinal direction^ They are represented by their projections 
in the x-y and x-z planes j hence, two planar curves are used to represent each 
three-dimensional curve. The parametric method of cubic splines (ref, 9) is 
used to curve-fit each planar curve, with the chordal distance between the 
coordinate points as the parameter. The parametric spline allows infinite slopes 
whereas the regular spline will not. To facilitate the application of the splines , 
slopes in the longitudinal direction are calculated at each cross section from 
the parametric spline. Then for the region between two successive cross sections^ 
the y and z coordinates of a three-dimensional curve are represented by cubic 
polynomials in x with the coefficients determined by. the coordinates and slopes 
at the two cross sections. When the slope of a coordinate at a cross section is 
infinite, then that coordinate is represented by a conic section for the region 
between that cross section and the one following it. If a longicudinal curve 
should not. give the desired variation, slopes at selected cross sections may be 
modified by the user, - 

Consider now the longitudinal variation of a conic section. In each cross- 

sectional plane, eq, (18) will hold but the coefficients A,, A„, A^, A,, and A^, 

1 4 3* 4 5 

will vary with ,x. As mentioned previously, these coefficients are determined 



by 4 defining points (the two control points # an intermediate pointy and the 
slope point) c For each segment, the 5 equations used to determine the coefficients 
Aq (q“l»‘'‘’°t 5) are formed by applying eq» (18) to the 3 points on the cross- 
sectional curve — the two control points y^, and y^^, , the intermediate 

point y^^, — and the slopes at the ends of the segment using the slope point 

MM (fc 

y , z (see Figure 6), This procedure yields the following 5 equations s 
s s 


where 


and 


pq 


5 

I 

q«l 


^ G A 

ill pq q 




„ _ 1 0 6 « C 


d = 
P 


-1 

-1 

-1 

-2 

-2 


-2 

— 

-2 

- 


^0 

y 2 

■'o 0 

z 

0 

^o 

z 

o 

-2 

^1 

hh 

-2 

"l 

^1 


-2 

^h 

Vh 

-2 

^h 

“h 



2z z 
0 s 



^^1^8 






(19) 


( 20 ) 


( 21 ) 


*The two equations for the slopes were combined with the first two equations 
(p = 1 and 2) to obtain the last two equations (p = 4 and 5)i> 
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At any longitudinal position j eqc (19) can be solved by any standard matrix 

inversion routine, such as SIMEQ on the CBC computer and SIMQ on the IBM computeY:;, 

- 2 

to determine the> coefficients A o • The derivatives dA /dx and d A /dx can be 

q ■ q q' 

obtained by differentiating eq® (19) and successively solving the resulting system of 
linear equations ^ The elements of G and their . derivatives with respect to x 

pq 

are obtained from the three-dimensional curves which Were spline-fit through the 
4 points used to define the conic section for that segment in each cross-sectional 
planer. 

Geometry in Polar Coordinates 

For some applications it is convenient to express the geometry inpolar 
coordinates x, f, and (j» where f = f(x,(^) Is the radius in a cross section measured 
from a reference point y^, , ieSo 

The ref erence ; point may be taken as y = 0 , z = 0 , but in some cases f is not 

IT 

single valued unless a reference point is chosen off the x axiso The angle i(i is 
given by (see Figure 8) 

<f> = tan~^[(z-Zj.)/(y~y^)] » 0 < ^ < 2ir (23) 

with ij) = ,0 corresponding to the line z “ z^ and y-y^ > 0 » Since y-y^ ~ f cos<ji 

and z-z^ = f sin<{> , eq^ (18) in polar coordinates becomes 

Tf^ + Bf + H = 0 (24) 

where .. 

2 2 

T cos (J) -h A2 cos((> sin(|) + A^ sin (f (25) 

B = (2A^y^ + 2A2Z^ + A^)cos((i + (2A2yj, + 2A^z^ + A^)sinq> (26) 

H = A^y^ + A^y^r + A^z^ + + A^z^, + 1 (27) 


15 



Figure 8 illustrates the polar coordinates in a cross sectional plane„ 

When solving eq« (24) for f from the quadratic formula, the + or ~ sign 
must be determined beforehand » To determine it, apply eq, (24) to the control 
point at the end of the segment (j-1), where f = f^ and (^ = (fi^ are known, to obtain 

-Bjf. ( 28 ) 

Substitute H from this equation into eq<> (24) and solve the resultant for f from the 

quadratic formula to get 1^2 

-B+[B^+4T(T fj + B f )3 ' 

f ^ J . , for T 0 (29; 

Now when and (29) gives ii tlifi sign is used when Bj+2Tjfj> 0 

and the - sign when B^+2T^fj <0 » At this longitudinal location, this same sign 
is applicable for any other value of <f> within the circumferential boundaries of 
this segment. In general, the sign must be determined at each longitudinal position.) 
The derivatives 9f/9x , 9f/9<j>, 92f/3x2 , 92 f/ 33 S^j, even higher deriva- 

tives, are obtained by differentiating eq, (24). Note, however, that the coefficients 
A 2 » A^, A^, and A^ are functions of x , These coefficients and their derivatives 
are determitied by the method described in the previous section, 

COMPUTATIONAL AJ.GORITHM 

Given the set of data points (y^,Zj^) in cross sectional planes at several long- 
itudinal stations, 

(1) For each cross-sectional plane, divide the data points into segments so 
that a conic section can be eurve-fit to the data points in each segment 
by the least-squares technique developed herein, 

(2) If the curves fit to the cross sectional data points are not satisfactory, 
modify them by one or more of the following methods; (a) define new 
boundaries (control points) for segments, (b) specify slope(s) at control 
point (s) (slopes may be finite or infinite and continuous or discontinuous), 
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(c) specify selected segments as straight lines, (d) a specific conic 
section can be specified for a segment 'by prescribing the slopes ac the 
ends of the segment and using only one. datum point between the end control 
points 0 , 

p) Represent the conic, section for each segment in a cross-rsectional plane in 
terms of ^ the two control points at the ends of the segment an interinediate 
point, and the slope point (see Figure 6)o 

(4) For each point found in step (3), spline-fit a three-dimensional lorsgitudinal 
curve through it and the corresponding points in other cross-sectional 
planes (see Figure 7)o 

(5) If the longitudinal curves are not satisfactory, modify them by one or more 
of the following methods $ (a) specify 8lope(s) at longitudinal station(s) 
(slopes may be finite or infinite and continuous or discontinuous), (b) 
specify selected longitudinal segments as straight lines, (c) redefine 

Che boundaries (control points) of the segments in the cross-sectional planes 
so that the points used for the longitudinal spline-fit form a smooth curveo 

(6) The geometrical . properties of the surface may be computed in polar coordi- 
nates at any position x, (|) by the following steps s 

a) Locate x between two consecutive longitudinal stations, and then 
locate the cross-sectional segment which contains 

b) Use the spline function to, calculate the coordinates, slopes, and 
second derivatiyes of the 4 longitudinal curves for this segment at 
X (see Figure 7 ) 0 

c) Calculate the coefficients Ap(p«l,^u 0,5) of the conic section at this 
location by use of eqo(19)> Determine the first and second derivatives 
of A with respect, to x ,from the first and second derivatives of eq, ' 

hr 

(19). 
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d) Calculate the body radius f frata ©q, (24) j, and the derivatives 
9f/3X( 9^f/9x^j 9 /9 g9^f/9x ?;* p arid eren higher 

derivatives from derivatives of eq<, (24) o 
(7) After a satisfactory surface fit has been obtained^ the data whlc.h nust bo 
retained for a geometry subroutine .package are the coordinates and icngi-= 
tudiaal slopes of the longitudinal curves at . those, longitudinal stations 
where cross-sectional data points were givenc Then^ ‘the geometrical 
properties of the surface can be calculated at any position by the methcd 
outlined in Step (6) above c 

Appendix C describes a computer program written for the IBM 360/175 computer 
to perform this computational aigocithmo This appendix also describes the input r 
and output data for the programs 

APPLICATION TO 70* DELTA WING 

The surface fitting method developed herein is applied to the 70® slab delta 
wing shown in Figure 9* The computer program described in Appendix. C was used to 
calculate the results p which are presented in Appendix Do This example was chosen 
because it .illustrates many of the, options available to modify the longitudinal 
curves and. because the results can be compared with an, exact solutiono 

Cross-sectional data are input at the 3 longitudinal stations shown in Figure 9 
andp due to symmetry ^ only the first quadrant is usedo Two segments (three control 
points) are needed. to represent the cross section at x«i0« The first segment is a 
straight line and the second is one-fourth of an ellipse o , The least-square curve^-'flt 
teehnique represents the ellipse exactly by specifying & *ero slope at control point . 

an infinite slope at control point j®3p and one datum point between theue two 
control points n Although two segments must also be used for th© other two croee •3sc* 
tionSp only one segment is necessary to specify the circle at .xc0c65798 and th<» 
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ellipse at .x»l»Oa Therefo re,, the first two data, points;, which are also control pointSj 
are made coincideato Than the exact curves are calculated from the leas t“square curve - 
flt-by specifying a zero slope at control point J“2s an infinite slope at 'contrel 
point and one. datum point between these two control points » . The exact location 
of -this datum point is irrelevant .except that it .must. lie on the desired curveo < 

Recall'that each three-dimensional longitudinal curve is represented by its 
projections in the x-y/and.x-z planes,, niarefpre, 16 longitudinal planar curves are 
used for this example (8 for each cross-sectional- segment) o Modifications , are made 
to the, initial spline fits to 14 of these curves since their . exact shapes are readily 
obtained from Figure 9« The initial spline fit to the curves for segment 1 and 
fox segment 2 are correct » and thus no modifications are necessary for these two 
longitudinal curves » 

The geometrical properties f and its derivatives are listed in Appendix -D for 
the circumferential angles ♦ 0, t/6 , x/3 , and ff/l at ^ ,0c3j ,1,0s 2,0s and 5, 

The results at x<°l,0, 2,0s and 5,0 are exact (yithln the accuracy of single precision 

< 

oh the IBM 360/175 computer) « whereas some, inachuracies are noted at :x«*0,3. Here 
the cross section is a circle, thus the exact solution gives 3f/3<|»»0s S^f/Sp^-oO, 
3^f/3xS^ “0, and the derivatives 3f/3x and3^f/3j^ should be independent of <t» f 
The surface fit could be made nearly exact at all positions by modifying the slopes 
of the longitudinal curves y^ and z^ for the second cross-sectional segment at 
x=>0» 65798, 

CONCL10H1H6 REMARKS 

An algorithm has been developed for surface fitting three-dimensional bodies. 

Data. points in several crpss-^sectienal planes of a body ere required as input data. 

The data points, in each cross section are then divided. into segments, and a least- 
squares technique Is used to curve fit portlpns of general conic sections t;o the 
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data points in each segment « The longitudinal variation of. the cross sections is 
established fitting parametx'ic splines through the coordinates of points us&d to 
define the conic sections in the cross'-sectlpnal planes <, If the surface fit is 
unsatisfactory ( it can be modified by specifying slopes at the ends of segments or 
requiring selected segments to be straight lines » The surface geometry can bast be 
analyzed through the use of an interactive computer graphics environment o 

After a satisfactory surface fit to the body has been obtained ^ data cards can 
be punched and used with some of the subprograms to form a geometry subroutine 
package^ This package is convenient for use in other computer programs requiring 
a mathematical .model of the body geometryo It will calculate the body coordinates^ 
slopes e and second derivatives e Computational time is .shorty and the amount of 
storage required is relatively small. The calculated data could also be used as 
input data for Coon’s patching method (ref. 2) » if the user should desire to employ 
that method in place of the one developed here« 
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APPENDIX A 


Solution for Conic Section £q.uatipn 

Substitute eqs, > (2), (3), (7), and (8) Into (1) to obtain the conic 
section In, the, form 

+ Pj2 + Qj • 0 (Al) 

where 

Pj “ [(A/m)j + (A/n)jly + (A/m)j <A2) 

and 

Qj * Ajy(y^y^) (A3) 

The solution of eq. (Al) is 

z - (-Pj + ^j^^)/(2C^) for Cj f ,0 (A4) 

where the discriminate is given by 

■‘i ■ "i * 

In order to obtain real roots for z from eq? (A4), eq= (A5) must givd 
Rj >0 for 0 < y < y^e 


First, note that eq»(AS) 

gives Rj > 0 when.AjCj > 0» ; Also 

, Rj » 

(A/m)j yj > 0 

at y “ 0 . 

and 



P m 

j 

(A/n)j 7y to 

at y « y^ 



Thus the possibilities are illustrated below* 




The tnin imuan value of A,C. which still gives real roots cotrespoads to the 

j .3 

case illustrated below 



aad dRj/dy «= o» The solution of » 0 and dR^/dy ° Q , using eq» . (A5) s, yields 


the minimum value of A.C . as 

J J 


(A^C^)* « (A/m)j (A/n)j 


(A6) 


Thus eq. (A4) will yield real roots for z if 

A^Cj > (A/m)j (A/n)j 


(A7) 


The choice of the sign to be used in eq» (A4) is deterndned by requiring 
the equation to satisfy th& coordinate of the control points s « y ^ 0 




and z “ 0 requires the - sign If (A/n). > 0 and the + sign If (A/n) . < Oo 
These conditions, are 2 dl compatible because a conic section passing through 
the two control points will have n^ < 0 if m^ > 0 and n^ > 0 if < 0« The 
sign given by these conditions can be used In eq. (A4) for all values of y in 
the range Q y 1 » 
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APPENDIX B 


Least-Squares Solution for Conic Section Coefficients 


Eq= (14) cannot generally be satisfied at all data points because 
there would be more equations than unknowns. Accordingly* eq. (14) is applied 
at data point k and rewritten as 


a. , A. , + B. , A . + Y. , A. .- + y, C. » r. , 

j,k j-1 "^j,k j 'j,k j+1 'k j j,k 


(Bi) 


where a. , * 3. r. » aud v values of a , 3., and y. evaluated 

jjK J J J 

at y « y, , z = z, in segment j; and> r . is called the residual. The least 

squares solution of the overdetermined system of equations determines the 

" Ti _ s t'lb. 0 .la ii 0 > IF-.iw s.txiioq .foi/jic.o ov?! hiIj 

coefficients , C’. , A. ,N) which minimize the sum of the residuals 

• ! ' • .L i I'.r .pe ni: bssu.jSc nao soor ; .yaoro tu jiyvlg ni;l8 

squared (see ref . 8) . 

. , 'f :■ Y 0 BgrxBl B.d j 

Define K, as the data point number which corresponds 'to the first control 
J 

point in segment j . Square eq. (Bl) and sum over all the data points in 
continuous segments to obtain 


N S+1 

I I 

j=l k=K^ 




N 

I 


K 


i+1 


j=l k»K, 






(E2) 


The right side of eq, (B2) is minimized by the system of equations obtained by 
setting partial derivatives of eq. (B2) with respect to the independent 
coefficients equal to zero. The result of setting partial derivatives with 
respect to equal to zero gives the following system of equations? 


a -A „ 
m,l m-2 

+ 

a - 

m,2 

"h 

^m,3^m 

+ 

a /A 

m,4 mrl 

+ ^ cA_. 0 
m,5 nri-2 


a - 

m,6 m-1 


%,7^m 


m,8 m+1 


m " 2, ' • • ,N (B3) 
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The result of setting partial derivatives with respect to 
zero gives the , following additional system of equations s 


,Gjj, equal to 


b ,A , + b „A +b ,A^, + b .C 
mtl m~l m,2 m m»3 Ttrfl. m»4 m 


m « 1, » ,N 


The combined system of equations, i.e.,, eqs. (B3) and (B4), gives a system of 
(2N')-1) linear : equations for (2N-1) coef f icients o The paremetexs used in eqs., 


(B3) and (B4) are defined as follows; 


K 


‘m.l 


m 


k*K 


, m-l,k m-l,k 
m'^l . * 


K 


m 


K 




^m,2 


^ y. ^m-ljk ^m-l,k , ^m,k ^m,k 

k*K , ’ * k“K ’ 


K 


m 


2 *^nri-l 2 \rt-2 , 

k«K . k»K k«K .. 


(B5) 


(B6) 


(B7) 


K 


nri’l 


*W2 

*m,4 ^ 'I ' ^m,k/^m,k | '*m+'l,k ^nri-ljk 

^ m ^“Vl 


(B8) 


^ntf2 

^tt,5 * ^ ■ °*m+l,k ^nH-l,k 


(B9> 


a X » ][ 

k-K 


^k \-l,k 


m-l 


(BiO) 


•Wi J 
m , 


(BID 
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K 


m,8 


I 


iiri-2 


a 


k»K 


k nri“l*k 


. 


ntfl 


b 1 

m,l 

“ Vl.8 

(313) 

^m,2 

* %,7 

(BIA) 

m,3 


(Bi5) 




^ j 

m,A 

If 

(B16) 


Note that if slopes should be specified at selected . control points,, then 
the least-squates solution described i4»ove is applied to the segments inbetween 
two consecutive control points with specified slopes « The minimum nundier of 
data points (exclusive of the control points) that, can be used in the least-*, 
squares solution is one in the first segment and two in each of the following 

segments o Any conic section can be prescribed for a segment by specifying 
the slopes at the ends of the segment and one datum point on the curve between 
the two control points. In this case eq, (BA) for m=l is the only equation to 
solve, and the least-squares technique determines the conic section which satisfies 
the specified conditions, provided a conic section can be fit through the 
prescribed slopes and datum point. 

In the solution of the combined system of eqs. (B3) and (BA), the terms 
and k interpreted as zero. Hence, the coefficients 

and appear in the resulting system, and recall that A^^ = 1 

unless the conic section should require Aj^ =' 0, 
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^PENDIX C 

Description of Computer Program 

Two computer programs are described here. The first determines the body 
geometry and allows modifications to. be made to cross sectional and' logitudinel. 

• f 

curves until a satisfactory geometry is achieved. The second program uses data 
cards punched from the first program along with some of , the subroutines to form, 
a geometry subroutine package for use in other computer programs. Both ^ computer 
programs are written lipi PORTBM IV for the IB>1360/175 compute^r* but with a small 
number of ; changes they can be used on the CDG ; coioputers « 

Program for Determining Geometry 

The, following subprograms are called from the main program or its sub" 
programss GEOM, ; INGEO, -LSTSQ, FZ, SPtX, SPLINE, FGN, DERIV, LINEQ, and-SlMQ. 

Geometry Subroutine Package 

After the PrPgrai^ for Determining Geometry has been modified to yield a 
satisfactory geometry,' cards are punched by setting NF«1 la the input data 
(see Description .of Input) . Then these data cards along with subprograms GEOM, .. 

FCN, DERiy, LINEQ, and SIHQ form a geometry;, subroutine package which can be 
used In other computer programs for detexming the body radius and its derivatives 
at prescribed locations ic and (|) . The statemeutd necessary ; to use this geomtry 
subroutine package in .another computer program are given below under the listing 
of the Main Program for Geometry Subroutine Package. 

Program for Determining Geometry , 

This main program calls subroutine INGBO which, in turn, calculates all the 
parameters necessary to describe the geometry, and punches this information on eards 
if desired. The body radius P and its derivatives FX, FP, FXX, FPF, FXP. are 
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calculated at prescribed locations x and , and the geometry is analyzed 
to see if it is satisfactory. Modifications are made to the gecmetry It 

is satisf abtory a 

Main Program for Using Geometry Subroutine Pact^ge 

The main program of some other computer program which uses this geometry 
subroutine package should have the COMMON and READ statements shown in the program 
listing,, A CALL statement to calculate the body. radius and its derivatives at 
prescribed values of x and is also listed. 

Subroutine GE0M(XX. PH, TR, ZR, E, EX. EP, EXX, EPP, Kg) 

At a given location (x, $) and reference point (y^, 5^, this subroutine 
calculates the body radius and its, derivatives. The arguments for this. subroutine 
are s 

XX* PH X, <|> 

ZR y^, 

F, FX, FF f, 9f/3x , 3f/3^ 

FXX,. FPP, FXP . 3^f/3,x^ , 3^f/3<|>^, 3?£/3x3<|> 

Other program variables arei 

CP, SP cosifr , sin<J» 

T, B, H . parameters defined by ,eqs. (25) , (26) , and (27) 

TX, TP 3T/3x , 3T/3<{» 

TXX, TPP, TXP 3^T/3x^- , 3^T/3 (|i 2, a2T/3x3<|i 

(similar definitions for BX^ BP, etc,) 

This subroutine is called from the Main Program. 
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Subreutiae INGEO 

This subroutine rea4s part of the input dats and calls other subroutines 
which read the retaainder of the input data » When prescribed* , it also punches 
cards with the, data necessary for the Geomet^ Subroutine Package =. Program 
variables are described in Description of Input and Output* and subroutines - 
LSTSQ* SPLX* FGN," and DERIV* .This subroutine is called from the Main Program 
for Determining Geometry. 


Subroutine LSTSQ. (114 m) 

This subroutine reads input data pertaining to a cross section and applies 
the le,astrsquares technique described herein to determine the -conic sections fcr 
each segment of the cross section. The duam^ arguments for this subroutine area 
IT longitudinal station hunger 

IiIY type of nose specified (see Description of Input) 

The variables involved with the READ and WRITE statements are defined in the 


Descriptions of Input and Output . Other program variables are 


A(I),* B(I), C(I), D(I)* E(I), GN(I) 
YY(IT* J, K) 


See Description of Output 


AM Aj/m^ in eq. (9) 

AN 

AMN (Aj/m^)(Aj/n^) in eq. (9) 

ALj bet, GA » and Yj , in eq, (Bl) 

BM(I) * BZ(t) Calculated yi z corresponding to input data points 


CAC AjC^ 

G(I,J) coefficient matrix for the system of eqs. (B3 ) and (B4) plus 

the equation A(l) « 1. 
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A<I) 


initially constant vector for the right side of the system 
of equations GX=A» After- calling SIMQ the solution vector 
X is stored in A» 


XC(J)j YC(J) y and z coordinates » respectively > of control point J 
ZM, ZN slopes and n^ in local coordinates 

If inequality in eq, (9) is not satisfied, the message CAC LT AMN is printed 
along with the value of AMN and I, Then is replaced by the value obtained 
using the equality sign in eq» (9)» 

This subroutine is called from subroutine INGEO onlyo 


Function FZ(XX. I. JJ. XXB. YYB) 

For a given cross-sectional segment and local coordinate y, this function 
calculates z(=FZ) in local coordinates and then transforms the point (y,a) to 
global coordinates (y, z)» Eq» (A4) in Appendix A is used to calculate z- The 
arguments for this function are 


XX 

value of local coordinate y 

I 

segment 

identification no. 

JJ 

segment 

number 

XXB 

y 


YYB 

z 



This function is used in subroutine LSTSQ only# 

Subroutine SPLX(NY,NJ) 

This subroutine fits a cubic spline through longitudinal points to form 
longitudinal curves. The initial spline fits may be modified by specifying slopes 
or straight-line segments at selected longitudinal stations^ See Description of 
Input Data for variables in READ statements. The arguments for this subroutine 
ares 

NY nose shape (see Description of Ihput) 
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NJ 


nundier: of longitudinal curves to be modified 


Other program variables are 

Y3a(IbJ,K) left-handed _BlGpe 

YXR ( I , J»R) right-handed slope 


see Description of Output 


3QP(l)p YP(I) parametric slopes dy/dS and ds/dS where. S is the chordal 
dis tance between data points 

This subroutine ,1s called from subroutine IN6E0 only<> 


Subroutine SPblNE (H,DEjK^AMjEpG, FP) 

rT ffi ■ iii^ iin •rnnrrtTTM r^^-^*— i ^' i m ■ — ■ -tmi ^ i <■> i ^ mihwik.i 

For an array of dependent variables this subroutine solves the tridiagonal 
matrix system of algebraic equations to detemdne the coefficients and elopes 
for the cubic .spline function as described in ref o. 9o The dummy arguments for 
this subroutine are 

array representing difference- between J*1 and J independent 
variables 

DE(J) array of dependent . variables 

K number of points, 1 < J < K 

AM(J) array of coefficients for spline function ('second derivatives 

at each point). 

E»G slope at the- first- and end points^, respectively 

N flag tp. describe end conditions of spline fits 

N « -2 indicates :A»(2) = AM(1) and AM(K-l) « ,AM(K) 

are used in place' of speclfy-ing E and Gi 

N “ ,-l indicates AM-(l) “- 0 and AM(K) 0 are 

used in place 'of specifying E and G 

N «= 0 indicates E le specified but AM(K-l) = AM(K) 

is used in; place of specifying G 

N « 1 indicates both E and G , are specified 

This eubroutlne Is .called from subroutine SPUC; only s 
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^uj^ro^^ine FCN_(I ^ J j>XX) 

This subxoutine calculates the y and z coordinates or 4 
by using a cubic variation in ^ between longitudinal stattuns^, 
derivatives of these coordinates are also calculated© The cuef ri^lci'/, - • 
cubic variation are decermined by the coordinates and slopes at Ltid u.s/v 
stations which surrouna the longitudinal position XX© If the slope ^ , 
nai station just aft of XX is infinite^ then the conic section 

Ay^ = AAx + BAx^ 
is used for this segment, where 

= y - y(l^l) and Ax -= X x(I-l) 

The coefficients A and B are determined by requiring the curve lo pass rn 
coordinates at stations I and I-l and have the prescribed slopes at these 
The arguments for this subroutine are 

I first longitudinal station which exceeds XX 

J cross-sectional segment number 

XX longitudinal position x 

Other program variables are 

YXR(l 5 J^K)p YXL(Is,Jj»K) See Description of Output 

YY(Ij,J,K) See Description of Oucput 

YL(LjMj,N) parameters which are calculated at x - XX 

L - 1 for y and 2 for z 

M - 1 far y (z if L ^ 2) 

o 0 

M - 2 foi y^ (z^ if L = 2) 

M - 3 for y. (z, if L = 2) 

n n 

M - 4 for y (z if L = 2) ' (see Figute '/ ) 

3 S 

N 1 for cc o r dmat G , 2 for derivative or d,, » 
respect to x, and 3 for second deiivaf i vc 

This subroutine is called from subroutines INGEO and GEOM, 
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Subroutine DERIV 

This subroutine uses the ■perafl»ter8 .YL<LjM,N) calculated in subroutine PCN to 
.calculate the .coefficients A^(p*>'l, ««. , 5) and their longitudinal derivatives . from 
eq» (19)0 Program variables are 

AA(,I) initially dp as glv^ by eq. (20) 5 ..after calling SIMQ it 

becomes Ac 
P 

AX(I), AXX(I) dAp/dx j d^Ap/di^ 

G(I,J) (see eq. (21)) 

GS(I(J) temporary storage for G(ItJ) 


GX(1,J)^ GXXdiJ) dGpq/dx , d^Gpq/dx^ 

YL(L,H,Np) see description in subroutine FCN 

This subroutine is called from subroutines INGEO and GEOM4 If, the conic section 

should be a straight line (the determinant of G^„ is aero^ and the solution is, obtained 

pq 

from subroutine LINEQ. 


Subroutine LINEQ (G«AA) 

This .subroutine solves the two linear equations 
G(lj4) X(l) + Gd.5) X(2) - AA(1) 

G(2,4) X(l) + G(2j5) X(2) » AA(2) 

for X(l) and X(2) and then stores the solution in AAd) and AA(2). The parameters 
AA(3)( AAC4), and AA(5) are set .equal to aero. If the determinant of the coefficients 
of the two equations above should be zero* the message *'DEN0M IN LINEQ IS 0" is ■ 
printed. This subroutine is called from subroutine DERIV whenever the conic section 
is a straight line. 

Subroutine ,SDtQ (G.A.N^KgS) 

This is one of the scientific ■subroutine packages in the IBM library to solve a 
system of linear equations GX»A. For tha CDC computer ( this subroutine may be 



replaced by SIMEQo The arguments for this subroutine are 


G(1»J) coefficient matrix In system GX»Ao 

A(I) vector A in the system GX “ Ac On r^eturn to calling prograr.,- r-he 

solution vector X is stored in Ac 

N the maximum order of G as stated in dimension statement oi 

calling programs 

K the order of G5 1 < K < N 

KS flag to indicate a solution is obtained or noto KS^l indicates 

no solution 0 
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Appendix C (Continued) 


BescTiptlon of Input 

The input data for the computer program are described io the. following step® 
’which, are also in the proper -sequence. 

Step lo 

READ. (1,100) NAM 
100 FORMAT (20A3) 

D 

D 

READ (1,102) NY, NX, NCT, NJ, NP 
102 FORMAT (515) 

The read statements for this step appear in subroutine INGEO. 

NAM' name used for body designation 

NY 1 for blunt nose, 2 for sharp pointed nose, 3 for cut - off nose 

(i.e. the body radius is non - zero at the nose) 

NX number of longitudinal stations where cross sectional s data are to 

be input. The nose station is included in NX although cross sectional 
data are not required there for blunted and sharp-pointed noses. 

NCT number of control points specified for all of the cross sections. 

Some control points may be made coincident .if a cross section should 
require less than this nund>er of control points. 

NJ number of longitudinal curves to be modified by specifying slopes 

or straight lines , at delected longitudinal stations . Note, that , 
there are 8 longitudinal curves for each segment in a cross section 
since the 4 three-dimensional. curves shown in Figure 7 are repre- 
sented by their projections in the x - .y and x - z planes. Also, 
note that the longitudinal curve through control points is common to 
adjacent segments, and modifications made to it must be input for 
both segments. 
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HP if cards are to be punched with the data necessary to use a geometJCy 

subroutine in some other computer program^ NP“1 c If cards are not v;o 
be punched NP/ 

Step 2<. 

The 4 following statetaents ere listed in subroutine INGEQo 
DO 2 1=1, NX 
READ (ljl03) X(l) 

103 FORMAT (F10.5) 

2 CALL LSTSQ(1,'NY) 

The statements below appear in subroutine LSTSQ <, 

READ (1,101) N,NLTj NST, N2T 
READ (1,101) (NC(J), J“l, NCT) 

101 FORMAT (1615) 

IF (NLT.GTo.O) READ (1,101) (NL(IL) , IL“1, NLT) 

IF (N2T.GT.0) READ (1,100) (N2 (12), DX2 (12) ,DY2(I2) ,I2'^1, N2T) 

IF (NST.GToO) READ (1,100) (NS (IS), DXX(IS), DY(IS), IS«1, NST) 

100 FORMAT (15, 2F10»5) 

READ (1,102) (XB(K), YE(K), K«1,N) 

102 FORMAT (8F10.5) 

0 

o 

e 

RETURN 

The read statements in this step pertain to the Information necessary to esteb’- 
llsh the curves in each cross sectional plane. 

X(I) value of x.at longitudinal station I, where I^l^ NX. 

The Information listed below is not req^lired for the first longitudinal station 
(I»l) on blunted and sharp - pointed noses. 

N number of data points for this cross section.. This number mesy vary 
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NLT.’ 


from one ' cTOss eectipn to anotlier > 
ntqnber of oegmexits to be specified as = straight lines., 

NST number of control points where left ~ handed slopes are to be 
specified 

N2T number of control points where right *- handed slopes are to be 
specified 

Note: The left ** handed. slope at a control point -is defined as the slope « 

In global coordinates (y»z),of-the segment which ends at the control 
pointf whereas the right - handed slope is defined as the slope of the . 
segment which; begins at, the control points If the slope at a 
control .point Is continuous, the left ~ handed r slope is the same as 
the right f- handed slopes A left - banded slope may be specified at 
control point NS (ISX as di/dy » DI(lS)/m(IS) , and this value is 
also used as the right ^ handed elope there unless a right- 
handed slope ris specified or a straight - line segment begins there > 

A right ,r handed slope may be specified at ^control point :N2 (12) as . 
dz/dy ^-DY2(I2)/BX2(12) provided a left - handed slope has already 
bem specif led at this same control point » Biscontinuous slopes 
at the beginning or end of a straight - line segment axe specified 
by dz/dy " 2)Y(IS)/BXX(IS) at control point NS(IS) « Slopes are 
specified in global coordinataa (y,s) by a, numerator and; a denomi- 
nator so that zero and infinite slopes may be input., The first segment 
of ,a cross section (J>1) must be either -a straight line or have the 
slope specified there [by di/dy « BYC1S)/DYX(XS)] . Also, ^ the last 
segment (J " NCTrrl) must either end. with a straight line or have the 
slope specified at the last control point (J * ItCT) by di/dy 
BY(IS)/BXX(IS) 
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NC(J) dati point nuaiber to be dasignatad a control point o The first snd l.art'r 
data points must be control points » and also 1»!C(J4'1) > li?CCJ)o 
number, of a segment -to be designated as a straight lina» 

N2<12) control point number where a right - handed , slope is specified. 

(see note above.) 

DX2(I2) denominator of right " handed slope di/dy at control point N2(I2)r, 
DY2(I2) numerator of right - handed slope da/dy at control point N2(I2)r, 

KS(IS) control point mpnber where a left,” handed slope is specified. 

(see note above) . 

BXX(IS) denominator of left - handed slope dz/dy at control point NS (IS). 
DY(IS) numerator of left - handed slope dz/dy at control point NS(XS). 

KB(K) y coordinate of Input data point K. 

YB(K) i coordinate of Input data point K. 

If this coa^uter program is to. be used with the body geomet^ry 
expressed in polar coordinates » the first data point (K 1) must be 
on the line $ >= 0 and succeeding data points should have 4(K+l)>i(K) 
relative to. the reference point y^ »z^ (see figure 8) . 

The read Ejtatements in this step read data for all the cross sections before 
proceeding to, Step 3. 

Step 3.. 

The following read statements are listed In subroutine SPLX and represent modi” 

f Ications to the longitudinal curves ° if no longitudinal curves are to be 

modified from the initial spline fits, then NJ “ 0 should be input in Step 1 

\ 

and the read statements in thia step will bd by - passed. 
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100 

101 

17 

J 

K : 

MLT 

MST 

M2T 

Notes 

MIT- 

ML(IL) 


IF .<NJ,EQ,0) RETUBN 
DO 17. IJK = 1,. NJ 

READ (1AO0) J.K, MLT> MST i M2T, MIT 
FORMAT (1615) 

IF (MLTiGT.O) READ (1,100) (ML(IL), IL “ 1, MLT) 

FORMAT (4(I1Q,F10.5)) 

IF (M2T.GT.0) READ (1,101) (M2 (12), DY2(I2), 12 =»1, M2T) 

IF (MST.GT.O) READ (1,101) (MS(IS), DYX(IS), IS l^ MST) 

IT (MIT.GT.O) BEAD (1,100) (MI(IT) , IT ® 1, MIT) 

CONTimiE 

segment, number in a cross - sectional -plane where 1<J^ (NCT “ 1) 
longitudinal 1 curve to be modified in segment J; 1 for y^, 2 for yj, , 

3 for 4 for y^, 5 for 2 ^, 6 for 2 ;, 7 for 2 ^^, 8 for 
(see Figures 6 and 7) . 

number of longitudinal segments to be specified as straight lines.* 
number of longitudinal points where finite, left - handedslopes- 
(dy/dx or dz/dx) are specified. 

number of longltudined points where finite, right - handed. slopes 
are specified.: 

Comments in the note after the description of N2T in Step 2 apply 

here also, except MS (IS) replaces NS (IS) , M2 (12) replaces N2(I2), and 

the slopes are (dy/dx or dz/dx) “ DYX(IS) at MS (IS) and (dy/dx or dz/dx) 
® DY2(I2) at M2 ( 12)0 However, infinite solpes are input separately .here. 

number of longitudinal . stations where infinite slopes in the long!”" 

tudinal direction are specified. . 

longitudinal station number where the beginning of a straight line 
segment is specified. The straight line will terminate at .the next 
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M2(I2) 

DY2(I2) 
MS (IS) 

DYX(IS) 
MI (IT) 


longitudinal station « 

longitudinal station number ^where a right handed slope is 
specified (see note above) 

right “ handed slope (dy/dx or di/dx) at . M2 (12) » 
longitudinal station number where a left “ handed slope is 
specified (see note above). 

left handed slope (dy/dx or di/dx) at MS (IS), 
longitudinal station number where the right ~ handed slope is 
specified as infinite. 
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Appendix 'C (Continued) 


Bescrlptien of Output 

The data described in Step 1 of the Description of Input are printed 
using the satae nomenclature , For each cross-sectiohal plane j, the data described 
in Step 2 of the Description of Input are printed with the same nomenclature 
except TB and ZB are the y and 2 coordinates of the input data points u Additional 
output for each cress section ares 

A,B,C,D,E cbefflclents of conic section given by eqto (1) 

SGI) sign to be used In eq<> (A4) 


SL0, SLl slope dz/dy in local coordinates at the beginning and 

end of a segment, respectively » 

YB(CALC) ZB (CALC) y and i corresponding to input data points, but cal- 

culated from conic section equation, 

YY(I,J,K) Coordinates of four points used to define conic 

section .in cross'-ssctional segment J at longi- 
tudinal Station X(I)^ . _ ^ 

K “ 1 for yo, 2 for yi, 3 for y. , 4 for y^, 

5 for Zq, 6 for zj, 7 for Zj^, 8”for z^ (see Figure 6) 


The data described above are printed for each cross sectional plane , Then 
the longitudinal slopes computed at each cross section by the method of splines 
(before any modifications are made), are printed, where 
1 longitudinal station number 

j cross-sectional segment number 

K same as K described In YY(I,J,K) above 

YXL(I,J,K, ) left - handed longitudinal slope, (dy/dx or dz/dx) 

YXR(I,J,K) right - handed longltudinsl slope. 


If modifications are made to the initial spilne-flt for the longitudinal 
curves, the data in Step 3 of Description of Input are printed with the same 
nomenclature, and thdn the modified left - handed (YXL) and right - handed (YXR) 
slopes are printed. 
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Next^ the coefficients for eq„ (18) are printed for each longitudinal 
station 1 and cross-sectional segment J» 

In the example used here, the following date are calculated and printed 
for prescribed longitudinal (x) and circumferential positions (<()) relative to 
the reference point 

X, PHI X, <t> 

FjFXjFP 3f/3x , af/3<j) 

FXX,FPP,FXB , 8 2f/3^2 ^ 
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Appendix C (Continued) 


Program Listing 
from IBM 360/175 Computer 

Contents Page 

Main Program for Beterminlng Geometry a o a ^ ^ « a e c o o o o o 44 

Main Program for Using Geometry Subroutine Package o c o » c ,o c o o .0 44 

Subroutine INGEC e 0 ® o e o » 0 0 » « e o o,© 0 © o © 0 0 0 0 o'o'©.& 44 

Subroutine SPX^ « a o.a ©‘0 ©.o q ‘5 © • • © «,* a © 6 o c'o o 6 q'& 6 . 45 

Subroutine SPLINE 0 0 « 0 0 « e o’© ©.• © © ©.o o'& o © © o 48 

Subroutine GEOM a a e 0 0 0 0 ©.© © © © 9 © ©,© © o ©.© n'o 4 ^ 

Subroutine FCN © «.o ©'© ©,© 9 » o e o © o © ©.© © 51 

Subroutine BERXV o © o « e.© o'© a » ©'e © © 0 © « © « ©,o o'o o o o' 51 

Subroutine LXNEQ ee »© ©.© o‘o 9« « • 0*0 ©a© 0 ©©e,© 0* 53 

Subroutine LSXSQ ©q a*©©.© o©©,© o« © © ©'© ©'© © © ©‘o©,© o' 53 

Function FZ a ©.o ©'© ©© ©© © © ©© o,© ©© & ,© ©■© © 0 ©o' 58 
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«JOB NCS.HAg.B2630/DEJAftNETTE»TlME»309f*A6ES«40,LlNES-5A 
C MAIN PROGRAM FOR DETERMINING GEOMETRY 
C EXAMPLE - 70 0E6 DELTA MING 
NRITEC3«666I 
666 FORMAT! IHU 
CALL INGEO 
WRITEOtlOO) 

106 FORMAT! ///fl7X«lHX9 LAX t3HPHl«lAX« IMF* lSXt2HFXt/tI8X,2HFP,13X,3HFKX 
l,l4Xf 3HFPPf 13X,3HFXP> 

DO 2 I»lt4 
X».3 

1F(I.EQ.2)X»1. 

IF(UEQ.3IX»2. 

IF!I*iQ«4$Xs3. 

DO 2 J»l»4 

PHs( J-n«3,141S92/6 

CALL 6EOMIX«PHfO«tO«*F«FX»FPfFXX«FPPtFXPl 
2 WRlTE!3tl0nX,PHfPfFX,FPtFXXtFPP*FXP 
161 FORMAT! />SXf4E16.7«/f9X»4E16.7) 

STOP 

END 


C MAIN PROGRAM FDR USING GEOMETRY SUBROUTINE PACKAGE 
e example « 70 0E6 DELTA MING 

C0MM0N/C0M2/X!20)fYY!20f 10t81iYXR!20tlOie)tYXL!20tlO»8lvYL!2o4«3)« 
INXfNCT 

RiAD!l»301)NX9NCr»NR 
301 FORMAT! 319) 

REA0!lt300)!X!n»l«l*NX) 

300 FORMAT! SE 16.7) 

REA0!l«300)!!!YYn«JfK)»Xal«8)t J>lfNR)tl»l»NX) 
RiAD(l«300)!!!YXL!I*J»K),K»1.8),J«l,NR)fI«2»NXI 

N|aNX«l 

RiAO!lt30O)!(!YXR!1«J«K)«K-lt8)fJ«l*NR)»l«lfNZl 
WRITE! 3 « 100) 

160 FORMAT! ///»17X>lHXfl4X93HPHI«14XflHF»19X92HFX*/«18Xt2HFP>l3X»3HFXX 
ltl4X«3HFPP*13X|3HFXPI 

60 2 I»L|4 
XXa.3 

tF<t»EQ«2)XXal, 

IF!l«i6«3)XX»2. 

lF!I«iQ«4)XXaS. 

00 2 Jel»4 
PHs01)*3« 141592/6 

CALL 6E0M!XfPH»0«*0«»F»FXtFPfFXX»FPPtFXPI 
2 WRITKBtlODXXfPHyFfFXfFPfFXXtFPPtFXP 

161 FORMAT! /f6Xf4G16.7t/»9K«4fil6.T) 

RETURN 

END ! 

SUBROUTINE 1N6E0 

e READ INPUT DATA AND CALCULATE COEFFICIENTS 

6 NY« 1 FOR BLUNT NOSE* 2 FOR SHARP POINTED NOSS« 3 FOR CUT OFF NOSE 
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C NX«NO« OF LONGITUDINAL STATlONSf NCT IS NO« OF CONTROL PT$.« NJ» NO« 

C OF LONGITUDINAL LINES TO BE MODIFIED 
C NP»1 FOR PUNCHING DATA CAROS 

CONN0N/CON2/X<2O)9YY<2O,t0»8l«YXR(2Ocl09B)>YXL(2O«10«82«YLI2ff4;3} ; 
iNXf^CT 

C0MM0N/C0M3/AI5I * AXI 5} »AXX( SI 
DIMENSION NAM(20) 

REAOIiflOOINAM 

100 F0RHAT<20A3) 

NRlTE(3tlOUNAH 

101 FORMATI 4X,20A3> 

REAOI If 102)NYtNXfNCTfNJ«NP 
102 FORMATISIS) 

WRITE(3f201) NYfNXfNCrtNJfNP 

201 FORMAT! /tSXf 3HNY«tt5t3Xf3HNX««I5f3X9AHNCT»fI5«3Xf3HNJa»15«3X93HN 
lP«fIS) 

DO 2 l»lfNX 
READIlf 103)Xm 
103 FORMAT! FID. S) 

2 CALL LSTSQllfNY) 

CALL SPLXINYfNJ) 

WRITE!3f202l 

202 F0RMAT!/,llXflHIt3XtlHJ»SXtAHAIl)*12X»4HA(2)fl2Xf4HA!3),/«AlX9 
1 4HAf4ltl2Xt4HAI5ll 

IB»1 

IF!NYfLE.2)lB«2 

NR»NCT-l 

00 3 I«1B»NX 

DO 3 J^ltNR 

CALL FCNUfJtXnn 

IF!yL(ltlfl)*EQ.YLIlf2»U.AN0,YL!2»lfU«EQ.YL<2»2«lMG0 TO 3 
CALL DERIV 

WRlTE!3t2031ItJtAmtA(2l9AI3l9A!4ltA15l 

3 CONTINUE 

203 F0RMAT(/t8X,2I4f3E16«7«/933X»2E16«TI 
IF (NP.NE.l) GO TO 4 
WRITE!2f301INXtNCT9NR 

301 FORMATISIS) 

HRITE!2«300l!X!I)fI»lfNX) 

300 FORMAT! 5E16.7) 

WRITE!2«300)(!!YY! IfJfK)9K»lt6) tJ«ltNR) fI»ltNX) 

WRlTE!2t300)!!!YXL! IfJfK)fK«l«B)9JalfNR)«Ia2»NX) 

NZ«NX-1 

WRITEI2»300)U!YXR! lyJtKltK-lfBlfJalfNRIfl-lfNZ) 

4 RETURN 
END 

SUBROUTINE SPLXlNYtNJ) 

C0MM0N/C0H2/X ( 20 ) t Y Y 1 20, IO 9 81 » YXR 1 20» 10, 8 ) t V XL ! 20 1 10 « 8 ) * YL t 2 « 4« 3 ) « 
IN, NCT 

DIMENSION YI20),YPI20|,XP(20),OS120),YM<20),XM!20),MLC5),M2!S)»OY2 
ll5),MS(5)rDYX<5),Mll5),XXI20| 

NR«N-1 

NQ-NCT-1 
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DO 2 
DO 2 

YUiaWOoJoKI 

DO 3 S-ggM ■ . ■ . . ' • ■ 

VUiaVVS . • . ‘ ‘ 

fiaKUH^U“l! . . 

T2“Vf 

3 DSn®£3SQ^?iTI=l‘«‘2472^*5«2} 

IF«f^Vo©Eo 2 IGQ TO 4 ■ ' • 

YPm = 0 

. SF«YY|2oJeKEoMEoQjypaiaOS<2im<2tJ9K» 

KP|U'-^0 

CALL SPLINEIDS?YeM 0 YMpWSllolo«O 9 YP) 

CALL SPLIi^ESDSv^»M«^»«KPnnLo»D»%P$ 

YKRUoJ9ttSsla,E20 
EF^YYSSoJsK^oEOoOIVXRa* JffKP^O 
IB =^2 

CD ?0 31 . • • . 

4 CALL SPLSNSIOS^YeNfYM,!* sl.^^atYP) 

CALL SPLWECDStXsNsKHsloslot-JjXPI 

• IB=1 
31 00 2 

VXLaoJ®K 5 =YPU J/XPU 8 
2 YXRUeJp!UaYXLU9J9KI 
WRS¥E092©3^ 

203 FOR^^AT1//0UX,lHI9 3XaHJ»3X9lHK»5XtlOHYXL^»JtK)96X»lOHYXR^»J9K» 

iSXslS^FOR I^IITIAL SPLINE! 

DO 32 J'^lsNQ 
WRSTEO940U 
401 FORMAT! 5X1 
DO 32 KSI 98 
I!RITE!39 558S 
558 F0R?^A?«§X1 

DO 32 I^ltW 

32 !i!RITE092@4l I 9 J 9 X 9 YKLI I » J 9 KI oYXR« I 9 J»K) 

204 F0RHAT!@X93I4»2E16«71 
IFfNJoEOcOmeTURN 

DO 17 XJK»l*N4 
TCUI^D 
MS«U®0 
M2UI^0 

MIC 0=30 

REAOU9l00IJgKcMLT9MST*M2T»MIT 

100 FORMAT! 1^155 

WRITEO 9205 I JsKsMLTsMST^MaTtNIT 

205 FORMAT!// 98 X 9 2HJ=, I393X*2HK=,I3,3X*4HMLT=,13f3X,4HMST®» I3t 3X* 

I 4m2Ti=(?|3j3X9 4HMIT=i.I3» 

IFIMLToGT<.0meAOU9 100MMLULl*lL«ltMLT) 

101 FORflATS4Sll09FlO«5n 

IF«M2ToGToQSREA0U9l01HM2U21f0Y2f 1:21 9l2»ltM2T) 
HFIMST®CTo0IReA0«l9 lOlHMSUSI 9 OYXIIS} ^IS^lrMSTl 
IFiMSf »@T<>0$R£AD(l«l00HHlfITlf ITslfNIT) 
IF!M2TcG7o0)NRI7EI39206MN2f I2)(0Y2n2)9l2ai»H2Tl 

.!F!^SnGT9OII!RITEI3e207HMSnS)9DYXnSlf IS«l»MSTl 
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tF(MIT.GT«0)WRITE(3*208HMIf ITItlT«l«MlT) 
IF(MlT.GT.0)WRlTE{3f209)(NL(lL)tlLeltHLT) 

206 FORMATi //« 10X,2HM2f7X,6HDY2/OX»/tC/f8X,I4«E16.7n 

207 FORMAT! //« 10Xt2HMS*7Xt5H0Y/0Xt/t l/»8Xtl4»E16«7) ) 

208 FORMAT!/ , 8X* SHMI® #514) 

209 FORMAT!/ «eX#3HML«f514) , 

1S**1 

: 12=1 

■ IT=l ■ „•/ . 

1B*1 

■ 00 5 I»l#MR ' 

OYY=YY! 1+lf J#K)-YY11»J»K) * - 

DXX=X!I+1)-X!1) 

OSS“SQRT( 0XX**2+DYY**2I 
IF!I.EQ.ML!IL))GO TO 6 
IFiUEQ.MKirnGO TO 7 
IF(I.EQ.M2!I2HG0 TO 8 
IF!I«EQ.MS!ISUGO TO 30 ^ ^ ^ 

IF!KA.EQ«nGO TO 9 

■ ' GO TO 5 :■ ■ .. 

9 MB=I 

V M8l=I+l ' *.' ■ 

KA=0 ■ 

DO 10 I2=MB1#N 
= ME=IZ 

.-r M=ME-MB+1 ; , . 

. . DYE=YY! IZ»J#K)-YY!IZ-l#JtlO 

OXE=X(tZ)-X(IZ-U 

DSE«SQRT!0YE**2+DXE**2) , ! 

IFlIZ.EQ.MSdSUGO TO 11 
tFnZ.EQ«ML!lL))GO TO 12 
10 CONTINUE : , V 

6 YXR!I#J.K)=DYY/DXX 
YXUI+1»J*KI=DYY/0XX 
IFd.GE.NRlGO TO 5 
DXE=X! I+2»-X!I+l) 

0Ye=YYd+2*J#K)-YY(m#JtX) 

DSE=SQRT!DX6**2*OYE*#2) 

XP! 1)»0SE/(DXE+0YY/0XX40YEI V>' 

YP!1>»XP! D40YY/0XX . 

YXR!I+1,J»K)=DYY/0XX 
y'- '■ KA*1 

lFUL«LT.HLT)n.= IL4>l 
lFd.EQ.MS!lS).AND«IS.LT.NST||S*>154l 
GO TO 5 

7 YP!U=DSS/0YY . - 

xp!i)=o - ^ 

YXR!I»J»KI=1.E20 ' 

IF(lT.LT.MITdT«lT+l • 

IF(1«EQ.MS11S)«AND.15.LT.NST)1$»1S41 

IB=2 

GO TO 9 

8 XP!l)=0SS/!DXX^0Y2d2)4DVY) 
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IF« I2ol?of-^2?n2“S2'»-i 
IF| iiSoLtof^STS SSofiS^a 

. GO T© f • . ‘ ' 

' S0 KPUa^0SS^3©^Jl<i‘®YStUS§^OW| 

VPSUaBVStd , - , • ' 

¥Xri St?JetU!=iVpf u/KiPS as 

SFnioL?of^sTSis=as<!*i 

GO TO ^ • ’ •' 

ll KP«SiJ|i^0SS/gPKE-&OVKMS8®0Yei 
'■ YP«M&C.0Y?I« 

. 60 TO as 

It PPB^SYYfi I24-i®JoKS“WU£5^ut^S»/4K«I£4-16“j:3I2D5 
XP3{4S£)0Sg/5 0KE-frFPS«^DVE J 

VPfiMSc3F^@OKP4Js^8 "■ • 

13 00 a^ a0^M®oi£ 

!R^--S@=n8<!>l 
YU^SoW IQo J©tU 

IPS a^oEQoUGQ TO 

TgsVYUQo Jg5^S-WOQ“a»JeK| 

OS? HKS3SflRTS?a00|*f200IS 

14 COWTIMUg 
• ■ IBai 

IF«WeiEo2SG0 TO m 

C4LL SPiSl^E?DSoYof«5s Yfle YP?U sYP? MS tleYPS 
CALI SPUft!E?0S5XXpPl9KM»KP«il9XP<MJ»l*XPJ 
l^ 00 IS 

YXL? JQo J(?{CS«YPO0S/XPaOJ 

15 YXR? J^e JoKS^^vpSiQS/kPUQl 
5 COWTIhSUi 

WRgfgOo^PIS JoKeJgS^ 

402 FORMAT? //oaiK9lHIe4X96HYXL«S9sa2elH#fl2fiHl95X96HYXRn®*I29 ^ 

ilFJsoIg© lH9e/S 

00 40f iQ-^loM 

403 WRSTEOsSiQI S^sYKiL? 109 J»KI ^YXfUIQ, J,K» 

210 FORMATS SXo I402il6o78 

IT COMTIMUi 
RETORf^ 
iMO 


SUBROUTSfiJg SPUME«Hs06eKjAM*E«6*N,FP) 

DIMENSION! ^9?HHOilR}»AM{K}«A<20)9B?20l»C(20l90(20l«Q(20UO(20}« 

lPS2®St?FPS2(DS 

C . E IS LEFT OER^gG SS RIGHT OiRVeNal FOR € 6 6 SPECIFIED* 

C N=0 FOR g SPEElFIgOgNs-l FOR NEITHER SPECIFIED AND M«0 AT BOTH ENOS 

C Ns“2 FOR M AT iN0 POINTS = M NEXT TO END POINTS 
ASIIsQo 

BSUasi„ 

CSISgioS ■ 


\ 
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Oai*3.*UDE(2l-DE( 1M/HI2I-EI/HI2I 
1F(N.UT,0)0( 1)«0. 

MKVa.5 

tFCN.E0.-UA(K»-0# 

lF(N.EQ.O)A(K)a*B(K> 

IF«N.EQ.-2)A(KI*-B<K) 

C(K)sO« 

0<K|«3.«( G-CDElK»-DECK-in/HIK»»/HtK> 
IF(N«LE.0)0IIO»0 

Qm«-cm 

U<1}«D(1) 

pm*BU) 

DO 20 J»2*K 
1F(J,EQ.K)G0 TO 21 
B(J)«1. 

CCJ)aH< J+U/CH( J)+H( J^m/2. 


A<J)=,5-CCJ) ' . ■ 

DtJ|a3,*<U0E(J4l»-0EIJII/HU+in-<DEIJ)-DECJ-UI/HCJn/(HU>+H(a^ 


iin 

21 P« JJaAt J)*Q< J-IHBI J! 
0(J»a-CCJ)/P(JJ 

20 U!J)*lDtJ>-A( J)*UIJ-1II/PU» 
AM|K)sU(Kl 
KM1»K-1 
DO 22 J»1«KM1 


JR=K-J 

AH( JR)=Q< JR>*AMC JR+IHUC JR) 

22 FP<JR+1)»<AM< JR+D/a+AHI JR)/6»*HUR+l)-MDE(JR+l)-0EUR)l/HUR+U 
FP<l)a<-AM(l)/3-AM^2l/6)*H*2KlDE«2)-DEI 1) )/H(2) 

RETURN 

END 


SUBROUTINE GEOM( XX*PH#VR#ZRtFtFX*FP*FXX,FPP,FXP) 

COMMON/COM2/X(20)f YY(20*1098)tVXR(20«10»8) tYXLC20« lO«8)«YL(2»493)e 
lIXtNCT 

C0MM0N/C0M3/AAI 5),AX(5)tAXXl5) 

C LOCATE I AND J FOR X AND PHI 
JX«NCT-l 

DO 2 I«2*IX ' 

IF(XX.LT.xn) )G0 TO 3 

2 CONTINUE 
I=IX 

3 DO 4 Jalf JX 
CALL FCNIItJtXX) 

0YJ»YLU,2,l)-YR 
OZJaYLl 2.2tl)-ZR 
FJaSQRT ( DYJ**2+OZJ**2) 

PHJaARCOSIDYJ/FJl ' 

IF<DZJ«LT»0)PHJ«2«0*3«141592-FHJ 

4 CONTINUE 
JaJX 

5 CALL OERIV 
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SPnOgJl/FJ • . . 

Bl« §o0OAAO5OVKogo©oAAfiIgCsIK-s-A&t4l 
B2o g o 0^=a A (! 2 P 0 W og o ®C>Ciia ^ i g AA 5 S I 

€Iq1£FO02 

Sgu|p>ooi , ■ 

S2^q2'C’SFCjCF 

TTqAAU SOe2«AAg SBOStFi^goMOI^Sg 

Ssi§j.oCiPOB2^SF - ... 

?3agio20?1‘‘^F<J 

!FdT3oSTo®SS©!(^a3. 

CF-eOSdFMD 

SPuSSi^gPH!) 

C2=^CP=vSOg 

S2oSP002 

Sgf^agosFOCF ■ . ' '■ ‘ 

C2F>3C2^=S2 

E&OCg-s-AAi! 2B*S2F'/2<!'AAnS!S*S2 

TKt^AKd li)CjCl-!>A?«{!21>OStF/2^AHiBJ«‘S2 ■' , ' 

lu^'Gs-i-AKKdgii^-^ssevg-^AKKjSKSg 
?P=f AAd 3B<’AAS U l)OS2F-3-Aa5 2S«^C2P 
TPPe^ 20 H & A (! 3 S - A A m HC 2I®“& A C 2 H S2P I 
T^P>o3a^OI<=AKd U BOS2P^AK!2J«^C2P 
BiK-^goOOAKd lS0YR<!-2o@0aK« 2?02R4&XKI 
S2Ka2o0OAXS2SOYK«2o0OAXS3?*$*IR-*-AX|5) 

BX!aBl^OCP<*B25tOSP> 

B I KK--^ 2o OOA KK« U om*2o QOA XX « g ) >S'Z R^A XX (41 
BgXX^goOOAXXdgBOm^SoOOAXKiSI^ZR+AXXCSJ ' . / 

BXX'--BIXXOCP*B2X^JOSP> 
epa-B|!&SFo®20CP 

SPpfi^iu 

BPXi^'-BlXOSP-s-BgKOe^ . ■ 

t2ciBOOg-<>o@OTOM 

TS‘^'0 

IFg?2o©To@i)TSaSm¥6T2B i 

IFdABSd T»oSToloi“@SB60 ?0 f • 

lP|BoMEo08F^°E/B 


Fn9-0-c-SGf»5OTSl>^2^T 

Oe-^or'O'iT'O-B 

P2:^F0iP 

PX^-id-PgOTX-FOBXB/Dg 

PSs^'.i{S^gopEsoogOf-4OFOppOT^-Pi0tPP“2<‘FpUi8P-F4iBI>P)/06 
PXX-- d -gOF X0020T-40P0F XOTX-Pg^f KX-24PX4BX-F4BXX ) /06 
yP'j^pgOfXP“2Op0FXOTP'-2OF«sFFW / 

20-roFpop^^pogPH-FK^.!BP“FFii^BX 
F8«>'-^gyFoyi^B^©y 
KtiTUPIW 
gs^p 



SU8RDUTINE FCN(I*J«XXt 
C CALCULATE X-VARIATION FROM SPLINES 

C YYn*J«K)« I AND J ARE GEOMETRIC POSITIONS K» 1 - YO^a - Ylt 
C 3 - YH» ^ - YS» 5 THRU 8 SAME FOR I 

C YL«LtM,Ml, i»l FOR Ye 2 FOR 2e M * I FOR YO* 2 FOR Vie 3 FOR VH, 4 FOR Y 
C N* I FOR F,2 FOR FX»3 FOR FXX 

C0MM0N/C0M2/X( 20 1 « YYI 20# 10#8t # YXRI 20«10#a) •YXLI 20# l0#8 1 #YL I2«4» 3} « 
IIX.NCT 

DIMENSION FI3*3) 

OXsXI n-X( I-l) 
u*(xx“Xii-in/Dx 
U2*U*»2 
U3*U**3 

FIl# 1)=2«U3-3<'U2*I 
F<2»U = (U3-2*U24U»*DX 
FI 3»II=(U3-U2I«0X 
FIl*2l»(6*U2-6PUI/0X 
FI 1# 3I«< l2*U-6) /0X*P2 
FI2,2I=I3*U2-4*U+U 
F(2#3)»t6*U-4)/0X 
FI3»2)»3*U2-2*U 

FI3»3»*(6*U-2»/0X / 

00 2 L»l,2 . 

DO 2 M=l#4 
K^A^IL-IKH 

lFtYXRII-lfJ*K).GE.l.E20)G0 TO 3 

DO 5 N» 1.3 

KK-1 

IFIN«GE.2)KK«0 

5 YL(L,M,N)*lYY(I-l#J.K)-YYIl»J,Kl>*FU#NKYYII,J#Rl*KK+YXRII-ltJ»KJ 
1*FI2,NI+YXLI 1#J.KI*FI3.N) 

60 TO 2 

C SEGMENT BEGINNING UlTH INFINITE SLOPE 

3 TUYYII,J*K»-YYCI-I.J.K) 

T*Tl/DX 

A=2,0*T1#I T-YXLI I. JfK») 

BaT*I2»0PYXLII. J,K)-TI 

SGN«l.O 

V=XX-XII-IJ 

IFITl«LT.0)SGN*-1.0 

DY=SGN*SQRT( A*V+B*V*#2 ) 

YLIL.M.UsYYII-l.J.KJ+OY 
IFIDY.NE.O) 60 TO 4 
YHL.M.2»«0*0 
YLILtM.3)B0.0 
GO TO 2 

4 YLtL,M.2>*IA+2*0*6*DX)/2.0/0Y 

, YL(L.M#3)»IB-VHLeM#2)4*2>/0V 

2 CONTINUE 

RETURN • 

END 

SUBROUTINE DERI V 

C CALCULATE AI I ), AI 21 #AI 3» . Af 4) # A|5» AND THEIR DERIVATIVES U.R.T.X 
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e0P-'iM0N(^'CDi^2/i'K(120U VVII20D IL0o@i) cVKRU 2©ei©(?8l gYKLS20ol0t)® I 0VLC2e4<?3J p 

• iPslpf^Cti' 

• pAK^dSS 

OEKIKMSSON QdSB pQKd SD o02UUgDp©0pS& o©KaSgSSo©Jl5U5o5J)o©S!foSS 
08 2 

QUB^3ViUSio4oE» 

@2^(!Jl^--VL(!JSp4o2 9 
JD-VlU3p4o3S) 

00 2 S^^iLoS 
©(!HaJI)=WJJ3pBpi8 
©HOo^D=^V0J J3pS p29 
©^^Hd SoJJaVILUSo Ho3&- 

DO 3 4“ 1. p 2 
■SPd JoS0o3SK=5 
EF 2 J oGE o 2 S L=S 
DO 3 E=Ep3 

©3 B p Jl B “©S S pK !) ®^G(i E p L S 

©JIS EpJ)0s6S EpKD=>gK2 IpLlloGJtd SoKS.&G« loLB 
f ©H22(l EpJS-©3 HyiU’^G^5?noLK206HOoK8»?‘©5«H9LS^6XXf Ip5C^«>&3 2pi2 
00 4 
E3aH“3 

©OpSS= 3G« ESoSS'SQS^S-O'SS I3p4D©Q2§5 

©ItOpSB^SUSpSDogKd^So-GHOSpSJ^Qd^UGI I3*4>*QX< 5)+6XU3c48>&gSSi 
T=3G3E3pSDC'0XUd45'O-2’^'©X8 E3p2H^QX2 4H0XX( I3«9)«Q(4) 

2p21)3lr^G(! I39^jl)OQHaSSD-!>2*GKi l3p4l«'QX|5>*GXXU3o4I^Q2 SJ 
DO 4 <J)t^4p5 
K - 1 

IFflJoE©c5i)K=.'3 

■ ©d SpKDr^SOGd E3pJ»0@UB ' ■ ■ 

©«S EpKl)^20d©«I39JD*QJ2d JSc-©X039J»*^QS J) » 

©KHS 2oKD= 32*S©2 E3oJ!D=5QH22MBo2C!®K5IBeJ84QXUHGXK039 J|!»0UU 
©d EpJlB-©(E HSpJHQdJi) 

©« fi Ip J) 2^^©Kd Oo B o© nd J D 

^ ©nKdSpJD^©K^3a3p^Doi9?^H(!J9 • , ■ 

©0 S 2^Ep3 

s ■ • . 

. DO & E-’Ep5 
©0 © J“Ep3 

€> ©S d K 0 sil B^Sd E 0 iJ 2 

C£i!L(L SlIflQ d ©^p AA 0 Sp ^ 9 2 

EFdKSog9oEi)©0 TT0 f ' . . 

BO T E^ Ep ^ 
ma E B-0 

00 1 ■ ' ■ ■ 

EJ'-AMd H &“AAdJSO(©|5dEoJ0 ' ' 

? ©SdSpJB=-^©d Hy^B 

CAf,E> S S N@ d ©Sp AK 9 !^o § pKS J 
E(Hd£tSof@oEO ©0 T© <S . 

©0 ® B^Epg 


\ 
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AKXU) = 0 
DO 8 J=l»5 

AXX(I)=AXX< n-AA( JJ^GXXf ItJ»-2*AXlJ|»GXUt Jl 
8 6 S< I«J)sG( 

CALL S1MQ16S«AXX«5«5«KS) ‘ 

iF<Ks<»eQoU GO TO 9 : . 

RETURN 

9 CONTINUE 
AA(U=»l 

AAC2)=-l ’ 

: CALL LlNEQtGeAAl 

AX( n»-AA(4)«GXn94l-AA( 5}*GXUt5) 
AXl2)s-AA<^I^GX«2«4l-AA(5J*6X<2*5> 

CALL LINEQCG^AX) 

AXXm=“AAC4J«GXXU94l-2*AXIAI*GXCl,4|-AAC5)«GXXa9 5J“2*AX<5J*GX 

m*5> 

AXX(2J*»AA«4J4GXXI2*4)-2*AXI4I*6XI2,4)-AA|5)*GXX<2*5>-2«AX«5»*GX 

l(2»5» 

CALL LINEQIG*AXXI 
RETURN 

- ■ END ^ • • 

' • 

SUBROUTINE LINEQIGfAAJ 
DIMENSION G(5,5)9AA<9) 

0EN»G(l*4)«G(2t5>-Gtl*5)*6<2»4) T 

lF<DEN«EQoO)GO TO 3 . 

UP4=AAU»*G(2»5)-GU*5J*AA«2I 

UP5«6U*4)*AA<2»“AA<U*6I2*4I 

AA(4»*UP4/DEN 

AAI5>»UP5/0EN ; 

DO 2 l = l»3 
2 AA( I|sO 
RETURN 

3 WRlTEi3,202) 

202 FORMAT! ///*2X,19H0EN0M IN LINED IS 0»//l 
RETURN 
END 

SUBROUTINE LSTSQ!IT,NY) 

C FOR THIS SUBROUTINE I IS USED FOR Xt X FOR Y* AND Y FOR Z 

C0MM0N/C0M2/Zt20lffYY(20»l0»S) tYXR! 20»10f 8 J « YXL (20» 10 » 8 I 9 YL C 2f 4» 3 ) j 
IIX*NCT 

DIMENSION NCC 10 S«XB( 20 JfY 8 ( 20 )»DXnO)fSDA( 10 >»CDAC 10 ) 9 CTC 10 ) 9 AGUO 
* l),B6U0)f ABC 10}«6G( 10) * BBUO) tAAC 11 ) » 6 EC 10 ) tBEC 10) 9 AEl 10 ) 9 EEC lOI^X 

ICZOIsYf 20)»GC20920)#BZ(20)*XSMO),YSClO)tXHC10)9YHUO)«BMC20) 
DIMENSION NU5) ,NSC 5)*0XXC5) tOYC5)*N2C5l,OX2C5)*DY2C5J 
C0MH0N/C0HI/XCU0)9YC( 10)»SAflO)fCAUO)»SGNC10)«A(20)9BUO)9CUO)« 
10(10)9EC10I ; . 

M=NCT“1 ; ' 

IFUT»EQ«loANDeNY«LE.2)G0 TO 32 V 

READII»101)N,NLT9NST#N2T 

REAOU* lODCNCC J)»J“ltNCT) 

101 F0RMATU6I5) 

. ■■■ NHI)«0 , . 
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M 2 S 1 

gP(!ML7o(ETo0mEA0II io 101 H Winns 

SF8Wl¥o6To0SKEfiiO(!lol©©5{!W2U2S@O)C2n2l ^OYSnSS p S 2 -loW 2 TS 
BP(IWS¥oGTo©S 1 ^ 1®0S fiWSfl 3SS ISS ISUISoloWSTS 

nm&^ It, u© 2 S gJJBgKS oVB«i;tJoS$!:>'5ioW3 
102 PO^KATg ©PlOo^S 

e W°¥07AL WOo OP D&TA IF>TSooWLn» SS WOo OF CONTROL PTo t?miCH IS BEGo OF 
C A -ifo LBWE SE©o9MiT IS TOTAL WQo OF OF ST. LINE SEGSo oS INILARILY 

C WSin mn WST ake for ewpot slopes «must have initial and end slopes 

C 0P^ HAVE A STo LEWS THERE! 

’ mETESSo^jO©! 

feiRKYSgSogm lum ' 
miTEOpaisD NoNlTdNsTpWst 
mSTggSo213HMC« EH I-lpNCT! 

SF 3 WLToCToOl)HRITE 8 3p21'?>! 5ML5 El o XnleHLTI 
IFgWSToGToOJWRITEg 3^2131 gWSU I pOKX{ II «0Yf I) 9 »NST I 
SFgW2Tc©To0n'dRETES3p216| gW2«D pDJt2U I «0Y2U I ,Iait,N2TI 
miTifiSpIlTH Eo»B« UpVBS E! 9 loi 5 ,ftsj 
III F 0 RKAT 3 / pSKpSMK =oE16o?S 

Z12 FORWATg/ p8Ko2HM=9O935to'^HWLT«0S393K84HNSTs9l353K«4HN2Ts0|3l 
21'4 PORWATg/ e©Mo3HWC=9l©I5S 
FORMAT!/ 0 8^9 3nWL^ 0 3E 5! 

21 s FORMAT! //olOHp 2 HWSo 8 i^o 2 HDV 0 ISJSpaHDEe/® « 8l^9i4«2E16o7H 

2 U’ FORMAT!//ol0Kol«9?J!o3H0¥29l3Ks3HD229/,I0X*I4»2E16.7H 

111 - FORMAT! //o any SHFToW©9?»o2HY®9i«e2HZB8/f I 8 X 9 I5« 2E l6o7l S. 
eo aa j==>i 9 we? 

JCaWCgJ! 

aa veuD.'-iVBst^D 
mHTiOo3®0l 

S@@ FORMAT! //9ll»olHJoiXo^HAUm2S?94HeCJ>tl2X,4HCUI»l2K94HD«Jl0/e2I.K 
ao^HiS^S9]l2Hoa-3S©W!.il!oil®!lo^MSL0Ul9lOX»6HSLlU.I I 
ES-1 

iL^a 

EtA-i 

12^1 

m U HQ-noM 2 

IF! S^oSQoWL! EU !©0 TO 8 

' SF! H0oSQoW2O2HOO TO If ' . 

SFU^oEGoMS! aSD 8©0 TO I® 

HF!stAoi0oa!)©o TO as ■ ■ .' 

©0 ¥0 11 
15 

KA->:j© 

©0 ® ss-n-mioWCT 
IFI I/oEQoWS! HSJ) ISO TO 
SFi IKoEQoWi! BL HOO TO M 
^ eowTBwyi 

i [3¥1”VS! HO^TB-Te! S^S 

m ne 3 Ho-c- a B- jc ! B© i) 

34 ■ ■ ' ' 

■ • ' • \ 



VYnT^IQ9l»=«CC IQ> 

YYUT» IQ*2> = XC(IQ+U 
YYUT,, IQ® 3l=t XCUQ)+XCUQ+in/2 
YYUT®SQ*4»=YY( IT,IQ»3I 
“ VYUT9lQ95» = YCnQl 

VY( ITj lQe6J=YCnQ^l> 

YY( IT? YCtIQ^lHYCUQ)l/2 

YY( lT®IQe8)=VYUT«IQ*T) 

IF( IL*LToNLTHL»JL+l 
IFC IL«tT«UIL=l 

IF( IQoEQeNSaSloAN0»!S«LT.NSmS»IS4'l 
KA«1 

AnQI=tO ' ' . 

B< IQI=0 
CUQ)«0 
01 1Q)=0 
E( IQ>=1 
< SON (IQ 1 = 1 
ZM=0 
ZN=0 

NRlTE<3®301ilQ®AnQ)tBnQI»CUQl90(IQ)9EnQ)«$GN(lQ}«ZN®£N 
K=NC( IQ I 

L=NC( IQ4-X)-1 ' . . 

0SB=S0RT(0Yl^«2-frDXl*«2J 

IF(DSB*LE«l.E“20IDSB»l 

CAA=DX1/0SB 

SAA=OYl/OSB 

00 33 JJ«K«L 

XHL = ( XB ( J J } - X6 ( K H AA« ( YB ( J J ) >y B ( K ) } ♦$ A A 
BM(JJ)»XC(IQ»^XHL«CAA ;• , •' 

33 BZUj» = YCUQ)-»-XHL«SAA 
GO TO 17 
10 OYl=OY( ISI 
DX1=DXX(1SI 
IF(lSoLT.NST5lS=lS4l 
GO TO 15 
19 0Yl*0Y2U2l 
DX1=DX2(I2J 
IF(!2,LT,N2TH2*I2+l 
IFlISoLT.NSTHS-IS+l ■ 

GO TO 15 ' ■ ■ ■ ; ■■■■ 

■ 12 ME=IZ ■ , • • ■. . 

OYN=YC« IZ+l»”YC( IZI 
DXN=XC« IZ-*-l>=XC( IZJ 

. M=ME“HB ^ • 

' 60 TO 7 , 

16 ME=IZ . -- .. 

0YN=DY«IS8 
DXN=OXX« ISI 
M=ME“MB 

' 7 00 14 I=lj20 
. ^ A(II=0 

00 14 J=*l®20 
14 G(le J|aO 
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V 


DO t 

TSi^ySSLD^-VBflKD 

OH3 Ti^o^oirioof i) , ■ ' 

eAos)=^Ti/0K(! n 
BFUoEQoaSGO TO 2 
SDA([ I! D-SAd EHe^CAS S«=12'-eAd 
COA« I j£3£Ad B 5^eAU=a KSAUS*SAO“ll 
eTOSsCOAfiES/SDAdO 

g COMTE wyg • • ' ‘ 

Ws d 0 Y !>>£ A S 1 E=D s A d 2. E J /d OK I A 8 U ■*• DV I *S A ! I H 

YPM=dDYN«'CAdKi“0KWOSAd KH/dDKW*CAI5idH'OYN*SA4MI i '‘ ' 
EF«Y^oe0o®oO^oTPAloE0o®^©0 TO ® 

CTUE"'--WP, 

CT(!M^EB“I^WM •, ^ ' 

SOAW»KS-a 

SOAdiB^a • ■ . ■ . , '■ 

DO 3 

ASdSEsO , ' 

0CdES-@ ' ' . ■ ■ . • ■ 

ADd a B "0 

C6 d E B S3 @ 

®0d E Bo© ' • 

AAdSB^O • ' ■ ' • - 

6EdIl)=5© . ■ . • , • 

©SdEB"© 

AEdE 3=0 ' ■ . ; ■ 

gSdS3=© 

K=MCd I^KB-E8 
L^MCd 

DO 3 JaKui 

KdJD®dKBd J 3“KB«ECJ5)^CAd'E3^dYBdJ3-Y8CK))«SAU| 

Ydy3^-dK®U3 = KBdKB J’J'SAd E 3*« YBU 3“YB4KH*CAt IJ 
ALsYUB^dDKdO^KdJBB/SOAdBB 

©Efs^Kd Jl^^^2-!>dCTd SS-CTd EoED gOKU&^'YUl-CTt I )«OXC n«'V4 43“©K02^KU3 

6A«Ug0S=^Yd JS/SDASI'{>18 

AC8 E 3^A6I 1 1 ■S’AL^GA 

06O3=BiBd ES-5-BET^GA 

AB(Un)=ABgS)oAL^0ET 

CSd S B=GGd I B '^SA^GA 

®©fi X 3=E0d E B*8-BE?^^2 •' •;.• 

AA d E 8 =AAd 1 3 ■• 

GgdE8=Gidn-3=GAe=Yd45!?«i2 ', . 

©Ed E8=D€d S J-s-SiT^SJYd 
Aid 1 8=A£3 1 3-S“Ai^5Y! J 
SEdI§=EEtIHY« 

3 COMTSMUE • ■ ■ 

AA3E^-S‘1J=@ 

Cdlo28=l 

Ada8=l 

Ilr-dKloi^ota©© ¥0 13 
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00 ^ !a 2 cH 

lFnoGEo3)GUol~2i«AGU-l) 

Gn 9 l-U=BG( I-1J+AB<I) 

GU^INGGU-lKBBlDMAf I + l) 
iFtSoLHoiM-lJ-lGf If l4U»BGirKABII>li 
!F« loVEolM-anGUsl+ZI^AGn+ll 
iMsUM:- 

GU 9 lM-ll*GE« 1 - 1 ) 

■ GnoIM)=BEU I 

IFUoLToMJGi It, IM-9-l)aAen + l) ' 

GUMt, l-U=Ae«U 
G< iM, n=BE( n 
IFn«LToM)GnM 9 l^l)= 6 Em 
4 GC IMjlM)=EE( n 
13 G«H^H 9 U=BEU) 

G«M+Iff2S=GE< U 

6 <H^l,M+l)=EEf U . 

I2M*2>9*M ■ ■■ ■. ■■ ' ■ . ' 

00 23 I = lpI 2 M \ ■ ■ ■ "V'V." 

DO 22 J=l® I2M 

IFfGUp J) oNEoOJGO TO 23 

22 CONTINUE / 

IFCHoEQoUGO TO 25 , , 

IZ=I+MB-1 

1F(5«GT«MHZ = U-M 
IFCNLT»LToIUGO TO 24 
DO 21 KI=ILpNLT 
2 i NUK^■s^UaMUKII 

24 NLT=NLT- 8 -l 

NL4 !Li»IZ . 

■ IF« IZoEQo!Q)GO TO 27 : ' 

GO TO 12 . ' . 

25 AIU = 1 ■ • - ' 

. A«2l*l/YP/YPN • V ' 

GO TO 26 

■ 27 NLtIL)»IQ^i ■ ■ 

GO TO 8 , ■ , ' 

23 CONTINUE . . 

CALL SIMQ!GpAt,20, I2 MjKS) if j: 

IF(KSoEQol) WRITE43,777) 

777 FORMAT* // 08 X 0 IBHNO SOLN FOR MATRIX) 

26 AM»1/YP ■ 

00 5 I=loM 

I: , '. 

IFU<.EQoUG0 TO 6 . 

AM=4A* I-l)~COASn*Am »/SOAU) 

6 AN=A(M)/VPN 

IFn»LToM)AN=*AU)=>CDAlI+l)-A(I+in/SOA(I+lJ 
CAC=A«M-fi-I)>>An ? 

AMN=AN*AM ■ 

IFICACoLToAMNSWRlTE<3*2l8) CAC*AMN#I 
210 FORMAT* ///pBXplOHCAC LT AMN*3X»4HCAC=*E16.7p4X«4HAMN»»El6«7p4X 
1 • 2HI»?i3p//l 
IF*CACoLToAMN!A«M4l | sAMN/A« II 
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I 


■ El 

Sf*3“AOI/AS^ 

IPI"AfSI/M 

KSie-gM^OKOB/fEM'i’EMl ■■•■.:■'-■ ' ;' ' 

‘ . ■ . • V ^ 

^SU l°KC^4l<‘XSi‘S'CAl IHYSi’^'SAOl . ' ■' 

VSU |sYCUHKSL^SaUl- 3 -¥SL«‘CAni ■■■ 

WLsFXgXSLeSsJsXHU'leVHUn ■ ' 

YYOTff Jt?U=XC« Jl ■.. 

VVUTo 492I = XCC J+ll 

■ YYUTovilpSI^^JCHUJ 
VYUToJff4l = XSn I 
WOToJ®5I^YCCJ» 

YYOfoJgSlssVCU^ll 
YYUf oJotlisVHm 
W«3 ToJ98I=»ysUI 

^^S?EO$30U JsA^ UdBUigCUHDUSfEUItSGNI lltZHoZM 
S@l raRKATI /t>8Xa24-t4il6«?j>^90Xtii6®?»F16.5#2El5«7) 

5C" !^£ 2 S “S’ P-l® i I 

00 5 JJ®KoL 

YHL==FIIKi JpBf?UJ®sBZUJ3l . 

3 COMTKE^UE 


0 I?HB£BVa?NI ' ' 

miTES® 92 l 9 H IpBMI SUBZUKI^lsNJ 

F0RHAIS//«8X«5HPT«N0t)5X»3HY5ICALC),BXt8HZB(CALC)«/t (8X»I5»2E16o7ll 
■ WI?E 0 3 ® 220 HT?ST?IT 9 ST 9 lT 9 ETpIT*IT 
21© F0Rf^&fC//9 UX,IH4,5K9 3HVVUE295H9 JjU*6X,3HYY<,I2,5H,J«2J,6Xj3HYYS 

IsIIsS^Jo J0 3lp6X9 3MYYl9l2eSHf Jsf4J s/*l8Xf 3HYYC »I2»5H» U»5U6Xe3HYYU 

I E205H9j06l06X93HYVC9!2v9H9J«7)»6X«3HYY(tI2«5HffJ«8>). 


2® yRITEOe22U 1»IYYI I?,| e«JU4«l»8l 

221 F0RSSfJST|/@gKf I4?4El&o79/el3X94E16«7) ■• . ■>. ' . • ,• '• ’ • . 

RgTURM ' ' , ... 

B 2 . 00 31 ’ 

00 3 1 . «J ^ ^ 0 8 ' ’ 

SI YVflToSgJSsO 
Riiua^ 

EMB ■ . '.. •. 

FUMCTIOW F£«XXf IsJJaXXBfYYSS 

€QmoS^/COKl/XCUOKYCUOUSAIIOI,CAaO)«SGNaO)«AI2Ol»BaOl9eU0l9 

101 10HE? 101 

TI-BI IHKX-^EUi ■ ’ • ^ ■ •' 

AU |«XX*!-D< I n 

EFIGi IKiGoOlGO 70 2 ' . ' ’ ' 

TS®W 2 ©Gm . ■ ’ 
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TA=T1#*2-T3 
IF(TA.LT*0»GO TO 3 
FZ={-Tl+SGN( n*SQRT(T<»n/CU)/2 
7 XXB»XC( JJ J+XX*CA(I )-FZ*SA< I ) 

YYB=YC( JJ »+XX*SA( n +FZ*CA( I) 

RETURN 
2 FZ^-TZ/Tl 
GO TO T 

3 URlTE(3e222l T4 

222 FORMAT! ///t 8 X 9 15W) I SCRIM INATE lSfE16.7«3X*lAHF0R FZ* SET *01 
FZ*-Tl/Cm/2 
GO TO 7 . 

END : . • 
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Input (Continued) 


I 4 1 0 0 1 

3 
1 

16 3 

1 2 3 

1 7 3 

1 2 3 

18 3 

1 2 3 

2 11 0 0 1 

3 - ; 

1 

2 3 1 0 0 1 

3 i 

2 4 1 0 0 I 

1 

2 5 3 

1 2 3 

2 6 2 0 0 1 

2 3 

1 

2 7 1 0 0 1 

3 
1 

2 8 2 0 0 1 

2 3 

1 ■' 



Appendix D (Continued) 


Outpujt 

70 06G DELTA WINS 

NY«b. 1 NX» 4 NCT= 3 NJ** 14 NP» 

X 0.6579800E 00 

N» 4 NUT* I NST* 2 N2T* 0 

NCw 124 


NL« 

1 


NS 

OY 

OZ 

2 

O.OOOOOOOE 00 

O.IOOOOOOE 01 

P 

O.IOOOOOOE 01 

O.OOOOOOOE 00 

PT.NO 

' YB 

ZB 

1 

0.9396926E 00 

O.OOOOOOOE OO 

2 

0.9396926E 00 

O.OOOOOOOE 00 

3 

0.6644630E 00 

0.6644630E 00 

4 

O.OOOOOOOE 00 

0.9396926E 00 

4 

A(J) 
E« JJ 

BCJ) 

SGNU) 

i 

O.OOOOOOOE 00 
O.IOOOOOOE 01 

O.OOOOOOOE 00 
1.00000 

2 

O.IOOOOOOE 01 
•0.1328925E 01 

oQ.OOOOOOOE 00 
-1.00000 

PT.NO 

YB«CALC) 

ZB(CALC) 

1 

0.9396926E 00 

O.OOOOOOOE 00 

2 

0.9396923E 00 

-0.3371769E-06 

3 

0.6644623E 00 

0.6644634E 00 

4 

O.OOOOOOOE 00 

0.9396926E 00 

J 

YY< 2fJ*U 
YY( 2*Jt5» 

YY( 2tJf2) 
YY( 2,J»6) 

1 

0.9396926E 00 
O.OOOOOOOE 00 

0.9396926E 00 
O.OOOOOOOE 00 

2 

0.9396926E 00 
O.OOOOOOOE 00 

O.OOOOOOOE 00 
0.9396926E 00 


C(JI 

SLO(J) 

O.OOOOOOOE 00 
O.OOOOOOOE 00 

0.9999935E 00 - 


YY( 2tJ^3) 
YYf 2,J,7) 

0.939692SE 00 

; o.oooooooE 00 

0.6644629E 00 
0.6644628E 00 


0(J> 

SLKJ i 

O.OOOOOOOE 00 
O.OOOOOOOE 00 

0.1328925E 01 
O.IOOOOOOE 01 


YY( 2tJf4» 
YYC 2,J»8) 

0.9396925E 00 
O.OOOOOOOE 00 

0.9396926E 00 
0.9396925E 00 
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X » O.IOOOOOOE 01 


N* 4 

NLT« 1 NST* 

2 N2t» 0 



NC» 

1 2 4 

• i ■ . 



NLa 

1 




O NS 

; DY 

02 



.'>■■■:, 2 1 

O.OOOOOOOE 00 

O.IOOOOOOE 01 

■ / ^ '• . ; 

'' 

3 

O.IOOOOOQE 01 

O.OOOOOOOE 00 



PT.NO 

YB 

ZB 


• 

1 

O.IOOOOOOE 01 

O.pOOOOOOE 00 



2 

O.IOOOOOOE 01 

0 . booooooE 00 

;■ ■■ . " 


3 

0.5000000E 00 

0.9216050E 00 



4 

O.OOOOOOOE 00 

0.1064178E 01 

f 



J 

A(J) 

/ 

B(J) 

Ct J) 

DCJ) 


E(J) 

SGNtJ) 

SLO(J) 

SLIU) 

1 

O.OOOOOOOE 00 
O.IOOOOOOE 01 

O.OOOOOOOE 00 
1.00000 

O.OOOOOOOE 00 
O.OOOOOOOE 00 

O.OOOOOOOE 00 
O.OOOOOOOE 00 

2 

O.IOOOOOOE 01 

0.1244837E 00 

0.1007751E 01 

-0.1460299E 01 


-0.1554019E 01 

-1.00000 

-0.9396920E 00 

0.1064178E 01 

PT.NO 

1 

. 2 
3 

.■■■' 4 

YBtCALCI 
O.IOOOOOOE 01 
0.9999997E 00 
0.4999997E 00 
O.OOOOOOOE 00 

ZB<CALC> 
O.OOOOOOOE 00 
-0.3240224E-06 
0.9216092E 00 
0.1064178E 01 




YY( 3fJ»l» 
YY{ 3»J*5J 

YY( 3»Jt2) 
YY( 3*Jf6I 

YY( 3»J#3» 
YYC 3,J,7) 

VYC 3,J,4) 
YY< 3,J,8) 

1 

O.IOOOOOOE 01 

O.IOOOOOOE 01 

O.IOOOOOOE 01 

O.IOOOOOOE 01 


O.OOOOOOOE 00 

O.OOOOOOOE 00 

O.OOOOOOOE 00 

O.OOOOOOOE 00 

2 

O.IOOOOOOE 01 
O.OOOOOOOE 00 

O.OOOOOOOE 00 
0.1064178E 01 

0.6886622E 00 
0.7716160E 00 

0.9999999E 00 
0.1064178E 01 


«t4t«s**«4<4i4f««;4c4i«4:)0e«4c««««4(«4t**4>4t«4c4>*4e4t^*4e)^*4‘*4‘4c«4t«4:4:« 


X« O.IOOOOOOE 02 

N» 4; NLT* I NST» 2 N2T« 0 
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OelOOOOOOg 01 
OelQQOOOOE 01 

oesooooooa 00 

OoOOOOOOOi 00 


Af J I 
EUl 

O^OOOODOOE 00 
OoiOOOOOOE 01 

O.IOOOOOOE 01 
"0«iS5402lE 01 


OoOOOOOOOE 00 
0o3275730E 01 
0o4l97335E 01 
0o433<9910E 01 


SCJJ 

O^OOOOOOQE 00 
lo 00000 

0«la44683g 00 

»lo 00000 


C( J) 

SLO( J) 

O.OOOOOOOE 00 
O.OOOOOOOE 00 

0.1007780E 01 
“O, 939691 IE 00 


ous 

SLIU^ 

OoOOOOOOOE 00 
0»OOOOOOOE 00 

»Oel460300E ©I 
0» 10641791 01 


YOfCALCS 
0«i000000E 01 
0s999999?E 00 
0e499999TE 00 
OoOOOOOOOE 00 


SB^CALCJ 
OoOOOOOOOE 00 
0«3275729E 01 
0«4197334E 01 
0o43399l0E 01 


t 


¥¥« 48J«U 
VV« 49J95> 

3«IOOOOOOE 01 
OoOOOOOOOE 00 


YYC 4«Js6l 

OclOOOOOOE 01 
0-32757301 Ql 


YV( 4,J,3) 
YY( 4tJf7l 

O-IOOOOOOE 01 
0-1637865E 01 


YY4 4«Js4J 
YYC 49J»88 

O-IOOOOOOE 01 
0-16378651 01 


OolOOOOOOE 01 O-OOOOOOOE 00 0.6866611E 00 O-IOOOOOOE ©1 

0.32757301 01 0-4339910E 01 0.4047346E 01 0.43399091 01 


J ^ ' YKtlI,JefC? YKRCr,J,K) FOR INITIAL SPUNE 


1 

I 

i 

1 


iUUUUUUUUUUUUUUUU 
1 Oc3i0756?E 00 
1 Qo 1068660E 00 

I “0.1074312E 00 


O.IOOOOOOE 21 
0-3187567E 00 
0.1068660E 00 
“0.1074312E 00 


I 


tuuuuuuuuuuuuuyyy odooooooE 21 



2 1 2 0«3187567E 00 
31 2 O 0 IO 6866 OE 00 
41 2 -0»1074312E 00 

11 3UUUUUUUUUUUUUUUU 
2 1 3 0«3187572E 00 
313 0«1068666E 00 
4 I 3 “0*1074317E 00 

1 I 4UUUUUUUUUUUUUUUU 
21 4 0«3187572E 00 
31 4 0»1068666E 00 
41 4 “0-1074317E 00 

1 / I 5UUUUUUUUUUUUUUUU 

2 : 1.5 O.OOOOOOOE 00 

3 1 5 OoOOOOOOOE 00 

4 I 5 0 »OOOOOOOE 00 

1 I 6UUUUUUUUUUUUUUUU 

2 I 6 -0ol954786E-02 
31 6 0o796426lE-02 
4 I 6 0«6901021E 00 

11 7UUUUUUUUUUUUUUUU 
2 17 -0ol070632E-02 
31 7 0-4355930E-02 
4 I 7 0«33l5780e 00 

I 1 8 UUUUUUUUUUUUUUUU 
21 8 -0*1070632E-02 
318 0-43559306-02 
41 8 0-3315780E 00 


12 tuuuuuuuuuuuuuuuu 

2:2 1 0o3187567E 00 

32 1 0-10666606 00 

42 I -0cl0743l2E 00 

I 2 2UUUUUUUUUUUUUUUU 

22 2 O-OOOOOOOE 00 

32 2 O-OOOOOOOE 00 

42 2 O-OOOOOOOE 00 

1 2 3UUUUUUUUUUUUUUUU 

22 3 0-1638253E 00 

32 3 0-2244857E-01 

4 2 3 -0-2192813E-01 

1 2 4UUUUUUUUUUUUUUUU 

2 2 4 0-3187567E 00 

3 2 4 0ol068656£ 00 

4 2 4 -0-1074307E 00 


0-3187567E 00 
0.1068660E 00 
-0-1074312E 00 

O-IOOOOOOE 21 
0.3187572E 00 
0-1068666E 00 . 
-0, 107431 7E 00 

O-IOOOOOOE 21 
0-3187572E 00 
0.L068666E 00 
-0.1074317E 00 

O-OOOOOOOE 00 
O.OOOOOOOE 00 
O-OOOOOOOE 00 
O-OOOOOOOE 00 

0. OOOOOOOE 00 
-0.1954786E-02 
0.79642616-02 
0.690L021E 00 

0.00000006 00 
-0.1070632E-02 
0.43559306-02 
0.3315780E 00 

O-OOOOOOOE 00 
-0.1070632E-02 
0.43559306-02 
0.3315780E 00 


O-IOOOOOOE 21 
0.3187567E 00 
0.1068660E 00 
-0.10743126 00 , ' 

O.OOOOOOOE 00 
O.OOOOOOOE 00 
O.OOOOOOOE 00 
0-00000006 00 

o.iooooooe 21 

0.1638253E 00 
0.2244857E-0i - 

-0.2192813E-01 

O-IOOOOOOE 21 
0.3187567E 00 
0.1068656E 00 ' 
-0.10743076 00 
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00 
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■ ®o795425l6- 

“®2 

©o796426IE»02 
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2 

$ 

0O&901O21E 
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00 
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■1 
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©oXOOOOOOE 

21 
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2 

6 

Qo4625672S 
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®o4625672E 

00 

$ 
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0o3a?3743E 

©0 

0o3l73745E 

00 


2 


0o4U7222E 

00 

0o4U7222E 

00 

a 

2 

Tuuuuuuuyyyuuyyyy 

Oolooooooe 

21 


2 . 

7 

0«3236437E 

0© 

0o3236457E 

00 

■3 

2 

7 

0S.3O887O3E 

00 

Oo3O087O3E 

00 

4 

2 

1 

0«4187793E 

00 

0o4l87793E 

00 

I 

2 • 

suuuuuuuuuuuuyuuy 

0<>1000000E 

21 

2 

2 

8 

0«4625650E 

©0 

0o462565OE 

00 

3 

2 

i 

©e3l73720E 

0© 

0o31r7372OE 

00 

4 

t 

i 

0o4U7240g 

00 

©o4U7249E 

00 

1 

■ Ka 

1 

HU^ 1 

S^STp 

0 H2Ts 

0 


S^Ijp I 
3 

I mu? U 1} VKRU? lo U 


g @«3639706E 00 

3 O^OOOOOOOE 00 

^ 00©OO0OOOE 00 


OoSasW^E 00 
0«0OOOO6OE 00 
-0al0743l2E m 


a K® 2 , 
E^ip a 

MLp 3 


I mi I? .1? 28 

ayyyoyuyuuuuuuuuy 

Z 0eB&B97D6g 00 

1 OeOOOOOOOg 00 

^ 0«000@000E 00 


I 0 ^2T= 0 


VKRUs i? 2 8 

@®IOOOOOOE 21 
@®36397©6g 00 
OoOOOOOOOi 00 
-®olOT4312E @0 


Jp I Ks S MLT® 1 f^STo 0 M2T« 0 . 
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MITa I 


MITe I 
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MI= 1 
ML» 3 


VKtUs le 3» 

• - -i ‘ 4'i 

lUUUUUUUUOUUUUUUU 
2 0e36397l4E 00 

'3 O.OOOOOOOE 00 

4 O.OOOOOOOE 00 


J« 1 K» 4 MLT= 
MI* 

ML* 3 


1 YXUI, 1, 4) 

I uuuuuuuuuuuuuuuu 

2 0.3639714E 00 

3 OeOOOOOOOE 00 

4 OEOOOOOOOE 00 


VKR(I» It 3t 

OolOOOOOOE 21 
0«3639714E 00 
OEOOOOOOOE 00 
“0 e1074317E 00 


1 MST» 0 M2T* 0 




YXRU, 1, 4) 

j 

O.IOOOOOOE 21' 
0«3639714E 00 
OeOOOOOOOE 00 
“0-1074317E 00 


J» I K« 6 MLT* 3 MST* 0 M2T» 0 
ML» 1 2 3 


1 YXL(I« 1» 6) YXRUt 1« 6) 

lUUUUUUUUUUUUUUUU OeOOOOOOOE 00 

2 O.OOOOOOOE 00 O.OOOOOOOE 00 

3 O.OOOOOOOE 00 0.3639700E 00 

4 0.3639700E 00 0.690102ZE 00 


J« 1 K* 7 MLTs 3 MST« 0 M2T* 0 
ML* i 2 3 


1 VXLUp 1s> 71 -- YXRU, 1, 7) 

j 

lUUUUUUUUUUUUUUUU O.OOOOOOOE 00 

2 O.OOOOOOOE 00 O.OOOOOOOE 00 

3 O.OOOOOOOE 00 0.181985QE 00 

4 0.18198S0E 00 0.3315780E 00 


MIT= I 


MIT*' 0 


1 

MIT* • 0 


V 
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j£3 I Ka @ 


3 © M2T° 0 MIT® 0. 


I 2 3 


■a .' mus i« 88 

puuuyuuuuuuuuuuu 
i oeooooooog 0© 
^ 9o@©@00OOE @6 
k @oS8|9850i 00 


Ja.. 2 I 

KI{? ’ 2. 

PlLci 3 


I mUe 29 IS 

'a yyyuuuuuuuuuuuyu 

£ yo3639706c 00 

i ©»©0©ooe©g 0© 

0«OOOOOOOE 00 


Jn. 2 Ka 3 MLfo . 

I / 

3 


E ‘ YXLUs ie 3S 

lyuuuuuuuuuuuuuuy 

2 0al726953E 00 

3 ®0ell92O93E'“O6 
“0oli92093B«>0ai 


Jo g 6 |e3 MLT~ 

I 
i 


E YKLUp 2^ 41 

a yyyuuuuyuuuuuuuu 

1 @*3fe39t0§i 00 ' 


radio Eo 88 



a ?dS¥o @ flgta © 


VKRd Ep' 2o 18 

@o 10000001 21 

0o3<&39706£ 00 
0o©©0000©g @0 
-©oi®?4312£ ©0 


I 0 M2T« • 0 


YKRUo 2o 38 

©olOOOOOOE 21 
Q«l?26953e 00 
-0»U92093E-06 
-O®2E9281BE“01 


1 KST« 9 nZT^ 0 


YKRUo 2s ^8 

0<»l©6O©0©g 21 
@®3fei0tO5S 00 


H17^ 1 


MITa 1 


HIT** I 



3 0.132454 7E-07 0. 1324548E-07 

4 0.1324548E-07 -0.1074307E 00 


J» 2 K® 5 MtT® 3 MST® 0 M2T® 0 MIT® 0 
ML* I 2 3 

> ■ } ■ ■ V ■ 


I YXUle 2« 51 YXR(I, 2, 5) A; 

lUUUUUUUUUUUUUUUU O.OOOOOOOE 00 ' 

> . 2 O.OOOOOOOE 00 O.OOOOOOOE 00 

3 O.OOOOOOOE 00 0.3639700E00 

4 0.3639700E 00 0.690l02ie 00 

'.-J® ' 2-. K« 6 MLT*' 2 , MST* 0 M2T®; 0-r ^ 1 ' 

• ■ ■ ' ^ 's i'V i : - 'V^? 

MI® , I . ' ■ 

■ HL® , 2 3 , /■■■,/■■■ 


1 YXL(l9 2» 6} YXRU* 2* 61 

lUUUUUUUUUUUUUUUU O.IOOOOOOE 21 

2 ,0.3639724E 00 0.3639724E 00 

3 0.3639724E 00 0.3639701E 00 

4 0.3639701E 00 0.4117222E 00 


J® 2 K® 7 MLT® 1 MST® 0 «2T* 0 MIT® 1 

'MI® I 
■ .ML® ■■ 3 


1 YXU lo 2s 7} YXRII* 2t 7J 

lUUUUUUUUUUUUUUUU O.IOOOOOOE 21 

2 0o3046873E 00 0.3046673E 00 

3 0.3639699E 00 0.3639700E 00 

4 0.3639700E 00 0.4I87793E 00 


J® 2 K® 0 MLT® 2 ' MST® 0 M2T® 0 MIT® 1 

MI® • , I . . ; ’ ■ 

ML® 2 3 
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. I ¥KLCl9 8 8 ■ VKROo 81 

luuuuuuuuuuuuuuuy ©oiaoooooe 21 


2 

3 

4 

0O 

00 

@® 

3639699E 00 
36397006 00 
3®>39t0li 00' 

0D 

@e 

©o 

36397O0E 00 
363970IE 00 
41172496 ©0 

■ 

I 


AliJ 


AC4) 

A(3) 

A{5) 

1 

1 

“0.11324T2E 

01 

-@o 95367436-06 
-Oo2©0674l6-O5 

-0.11324766 01 
0.32359746-06 

3 

I 

-O.lOOOOlOi 

01 

-®o8583069E“0S 
• ©o 97977356-05 

-0.88303206 00 
0.10887166-04 

4 

1 

©oOOOODOOE 

00 

OoOOQOOOOE 00 
-©.10000006 01 

O.OOOOOOOE 00 
O.OOOOOOOE 00 

4 

1 

0oll79822i 

0© 

0.26464466-04 

0.1041781E 00 


=@oU44409E-03 -0.6825432E 00 


0cS23S98?E 00 
“0ol754072E 01 

0®I047197E @1 
“0ol7S4CmE 01 


F 

FPP * 

0-7141427E 00 
0.6284285E-01 

0,7185342E 00 
“0.2382276E-01 

0«7185337E 00 
“0«238218lE-0l 

0*7141420E 00 
0,6284434E-01 


FX 

FXP 

O.4709243E 00 
-0.nil489E-0§ 

0.4492615E 00 
“0. 6265515E-01 

0.4492621E 00 . 
0«6265879E«-01 

0-4789269E 00 
0.3704963E-0S 


X 

FP 

OeBOOOOOOE 00 
*^0o4168O83H=O6 

0o3©OOOOO€ 00 
0b90152?OE~O2 

©«30OO0OOE 00 
-0o901661lE“02 

©oSOOOOOOE 00 
°@<,3241843E-06 

©oiOOOOOOE 01 
@oUS204§e°0S 

OolOOOOOOE 01 
©®S293887E“01 

©olOOOOOOE 01 
58133396-01 

OolOOOQOOE 01 
— 0 » 33143386— 04 

20000006 01 
@a32l4021i 0® 


PHI 

FXK 

O«OOOOO0Oi 00 
-0el721457E 01 


■O*15707946 01 
-0^17214556 01 

OoOOOOOOOE 00 
-0ail69780S 0© 

0^52359876 00 
-0«91?27i46-0l 

0®lO47l9?e 01 
-0e335627iE“0l 

0<,137©7966 01 
0e4245l736-04 

0.00000006 0© 
-©.14098O2E 00 


0.9999990E 00 
0.1169729E 00 

0.1014950E 01 
0.69441B5E-01 

0.1046981E 01 
-0.5744129E-01 


-0.2955627E-06 
0.3213953E 00 

0.16553786 00 
0.30399466 00 

0.30510246 00 
0.2100322E 00 


0.10641776 01 0.36397026 00 

-0.14098236 00 0.94669126-06 

0.93969166 00 -0.1244882 6 00 

0.21982266 00 0.32140306 0@ 



0« -Oft O* 


'' 2000000E 01 0.5235987E 00 

•4050109E 00 -0«1057327E 00 


2000000E 01 0«104?197E 01 

.3244193E 00 -0*3524287E-0l 


2000000E 01 
.3481045E-05 

5000000E 01 
.4086256E 07 


0.1570796E 01 
0«1960301E»06 

O.OOOOOOOE 00 
0.7555792E 14 


0oll33538E 01 
0.5557603E-01 

0.133S164E 01 
-0*4042364E 00 

0.1428X46E 01 
»0,7423788E 00 

0.4032016E-04 
0«1316266E 25 


5000000E 01 
•18270276 07 

5000000E 01 
•2170684E 01 


0.5235987E 00 
0.2014248E 18 

0* 10471976 01 
-0-9872645E-01 


0.6621871E 00 
0«7795724E 19 

0.19511166 01 
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V. 



LEGEND 


1. ’ < (A/m)i(A/n)j hyperbola which gives complex 

^ ^ ^ roots for z 

2, A,C. * (A/m).(A/n). ; _ v-:- 



Figure 4. Effect of A^Cj on conic section. ' 
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